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THE  LOT  SCHEDULING  PROBLEM 
IN  THE  HIERARCHY  OF  DECISION  MODELS 
Robert  Sheldon,  Ph.D.,  Capt,  USAF 
Cornell  University  1986 

The  economic  lot  scheduling  problem  (ELSP)  surfaces  from  competion  among 
products  for  a  scarce  resource,  usually  machine  time.  When  scheduling 
production  of  batches  in  this  environment,  two  issues  must  be  resolved:  the 
size  of  batches  and  the  start  times  for  production  of  each  batch. 

Standard  approaches  to  the  ELSP  focus  on  scheduling  multiple  products  on 
a  single  machine.  This  thesis  addresses  three  important  issues  that  put  the 
scheduling  problem  in  the  context  of  its  physical  setting  and  range  of 
parameters:  idle  time,  the  zero  switch  rule,  and  stochasitc  input  to  a 
bottleneck  machine. 

In  most  scheduling  heuristics,  the  reason  for  idle  time  is  to  balance 
cyclic  production  patterns.  Idle  time  is  also  optimal  in  solutions  to 
problems  with  high  setup  costs.  We  show  that  the  condition  for  inducing  idle 
time,  given  zero  setup  costs,  is  when  one  product  has  dominant  holding  costs 
and  the  remaining  products  have  low  machine  utilization. 

A  common  policy  in  scheduling  is  to  start  production  only  after  the 
inventory  reaches  zero.  This  policy  is  called  the  zero  switch  rule  (ZSR)  and 
is  regarded  as  a  good  scheduling  policy.  We  show  that  the  condition  when  ZSR 
is  not  optimal  is  when  the  ZSR  solution  yields  lumpy  production  patterns  for  a 
product  with  dominant  holding  costs. 

The  standard  approach  to  scheduling  considers  the  input  process  to  be 
deterministic  and  ignores  delivery  of  raw  parts.  This  thesis  examines 
scheduling  a  bottleneck  machine  with  stochastic  inputs  under  a  variety  of 
situations. 
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THE  LOT  SCHEDULING  PROBLEM 
IN  THE  HIERARCHY  OF  DECISION  MODELS 

Robert  Sheldon,  Ph.D. 

Cornell  University  1986 

The  economic  lot  scheduling  problem  (ELSP)  surfaces  from 
competition  among  products  for  a  scarce  resource,  usually  machine  time. 
When  scheduling  production  of  batches  in  this  environment,  two  issues 
must  be  resolved:  the  size  of  batches  and  the  start  times  for 
production  of  each  batch. 

Standard  approaches  to  the  ELSP  look  only  at  scheduling  multiple 
products  on  a  single  machine.  To  put  the  scheduling  problem  in  proper 
context,  we  examine  how  scheduling  decisions  on  a  machine  affect  and 
are  affected  by  other  decisions  and  the  physical  structure  of  the 
system.  This  thesis  addresses  three  important  issues  that  put  the 
scheduling  problem  in  the  context  of  its  physical  setting  and  range  of 
parameters:  idle  time,  the  zero  switch  rule,  and  stochastic  input  to  a 
bottleneck  machine. 

In  most  scheduling  heuristics,  the  reason  for  idle  time  is  to 
balance  the  cyclic  production  patterns.  Idle  time  is  also  optimal  in 
solutions  to  problems  with  high  setup  costs.  We  show  that  the  condition 
for  inducing  idle  time,  given  zero  setup  costs,  is  when  one  product  has 


dominant  holding  costs  and  the  remaining  products  have  lov  machine 
utilization. 

A  common  policy  in  scheduling  is  to  start  production  only  after 
the  inventory  reaches  zero.  This  policy  is  called  the  zero  switch  rule 
(ZSR)  and  is  regarded  as  a  good  scheduling  policy.  Ve  show  that  the 
condition  when  the  ZSR  is  not  optimal  is  when  the  ZSR  solution  yields 
lumpy  production  patterns  for  a  product  with  dominant  holding  costs. 

The  standard  approach  to  scheduling  considers  the  input  process  to 
be  deterministic  and  ignores  delivery  of  raw  parts.  This  thesis 
examines  scheduling  a  bottleneck  machine  with  stochastic  inputs  under  a 
variety  of  situations.  First,  we  isolate  the  issue  of  scheduling 
deliveries  to  a  machine.  We  look  at  using  state  information  to  schedule 
the  pre-bottleneck  machines.  Next,  we  consider  an  aggregate  planning 
model  to  evaluate  both  lot  scheduling  and  delivery  and  develop  an 
algorithm  for  solving  this  problem.  We  show  that  the  conditions  when  we 
should  consider  the  delivery  issue  in  conjunction  with  the  lot  sizing 
issue  are  when  the  holding  cost  of  raw  parts  is  high  and  when  the 
variance  of  delivery  times  is  high.  Then  we  examine  a  dynamic 
programming  formulation  and  consider  a  variation  of  the  assembly  model. 
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CHAPTER  1 


INTRODUCTION 

The  lot  scheduling  problem  surfaces  from  competition  among 
products  for  a  single  scarce  resource,  usually  machine  time.  Batch 
production  often  is  a  natural  consequence  of  manufacturing  several 
products  on  the  same  machine,  unless  setup  times  are  negligible.  When 
scheduling  production  of  batches  in  an  environment  having  a  single 
constrained  resource,  two  issues  must  be  resolved:  the  size  of  the 
production  batches  and  the  start  times  for  production  of  each  batch. 

The  solution  to  this  problem  is  computationally  difficult  (Hsu,  SO)  and 
therefore  practical  solutions  must  be  obtained  using  heuristic 
techniques.  This  scheduling  problem  occurs  so  frequently  that  it  has 
led  to  a  standardized  characterization  as  the  Economic  Lot  Scheduling 
Problem  (ELSP). 

Standard  approaches  to  the  ELSP  look  only  at  the  issue  of 
scheduling  multiple  products  on  a  single  machine.  However,  in  a  more 
realistic  setting,  we  must  ask  what  is  the  real  problem  we  are  trying 
to  solve?  Couched  in  this  setting  is  a  hierarchy  of  decisions  that  must 
be  made  at  different  levels  of  management  at  different  points  in  time 
and  at  different  locations  within  the  manufacturing  system.  Hence,  to 
put  the  scheduling  problem  in  the  proper  context,  we  must  examine  how 
the  scheduling  decisions  on  that  machine  affect  and  are  affected  by 
other  decisions  and  the  physical  structure  of  the  system. 
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The  hierarchy  of  decisions  for  most  manufacturing  environments  can 
be  broken  into  the  following  four  levels: 

Manufacturing  Systems  Planning  (or  Capacity  Planning) 

Production  Planning 
Plow  Planning 
and  Scheduling. 

The  decisions  made  at  these  various  levels  range  from  design  of 
facilities  down  to  real  time  detailed  scheduling.  Any  modelling 
approach  to  a  manufacturing  system  should  therefore  address  how  the 
model  fits  into  the  framework  or  structure  of  the  overall  decision 
process  as  well  as  how  the  problems,  issues,  and  decisions  from  that 
model  interact  with  the  rest  of  the  system. 

The  standard  approach  to  lot  scheduling  considers  the  input 
process  to  a  machine  to  be  deterministic.  Therefore,  the  issue  of 
delivery  of  these  raw  parts  is  ignored  when  scheduling  the  machine.  If 
the  delivery  of  the  parts  is  indeed  deterministic,  then  deliveries  can 
be  scheduled  to  arrive  just-in-time,  and  hence,  the  scheduling  of  a 
machine  can  be  looked  at  independently  of  its  predecessors.  However,  if 
the  inputs  to  a  machine  are  stochastic,  that  is,  processing  times  on 
the  predecessor  machines  are  variable  or  the  overall  delivery  times  are 
inconsistent,  then  the  schedule  of  the  machine  should  be  developed 
considering  the  issue  of  delivery  of  raw  parts  to  the  machine. 

A  traditional  ELSP  assumption  is  that  the  demand  process  is 
constant  and  continuous.  An  important  issue  in  the  physical  context  of 
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the  machine  environment  which  deviates  from  this  assumption  arises  if 
all  different  parts  produced  on  a  machine  are  subsequently  assembled. 

In  the  assembly  model,  the  work-in-process  (WIP)  inventory  is  held 
until  all  parts  are  available  for  assembly  in  addition  to  the  normal 
inventory  accumulated  due  to  batch  production.  In  this  case,  we  get  a 
different  inventory  pattern  from  the  traditional  ELSP  sawtooth 
inventory  pattern.  This  requires  a  different  approach  for  the  assembly 
model  from  the  tradtional  ELSP  approach.  Another  variation  occurs  if 
the  demand  process  is  dynamic  and  backlogging  of  demand  and  machine 
capacity  is  permitted.  These  problems  can  be  addressed  in  a  general  way 
for  arbitrary  ranges  of  parameters,  however  for  restricted  ranges  of 
the  parameters,  a  different  method  may  work  better. 

Traditional  approaches  to  the  ELSP  use  the  same  solution  technique 
for  all  scheduling  problems  regardless  of  the  range  of  the  parameters 
involved.  If  we  place  the  problem  in  the  context  of  the  parameters 
involved,  an  important  concept  to  look  at  before  addressing  any  given 
scheduling  problem  is  the  notion  of  a  dominant  product.  By  this  is 
meant  that  one  of  the  products  has  its  parameters  such  that  any 
solution  to  the  problem  will  always  be  dominated  by  that  product.  If 
this  is  the  case  for  a  given  problem,  we  can  focus  in  on  solution 
approaches  that  take  this  into  consideration.  Using  this  notion,  we  can 
get  better  solutions  for  the  restricted  class  of  problems  without  the 
added  effort  of  a  generalized  solution  technique. 
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The  concept  of  a  dominant  product  plays  heavily  in  looking  at 
whether  to  induce  idle  time  into  the  lot  scheduling  problem.  In  most 
heuristics,  the  reason  for  idle  time  is  to  balance  the  cyclic 
production  pattern  to  accomodate  non-rotation  cycles.  Idle  time  is  also 
optimal  in  solutions  to  problems  with  artificially  high  setup  costs, 
that  is,  the  cost  of  setup  is  higher  than  the  imputed  value  of  lost 
machine  capacity  due  to  the  setup  time.  Other  than  these  cases,  if  we 
assume  zero  setup  costs,  then  an  optimal  schedule  would  rarely  have 
idle  time  unless  one  product  is  dominant. 

A  common  policy  in  scheduling  is  to  start  production  of  a 
particular  product  only  after  the  inventory  of  that  product  reaches 
zero.  This  policy  is  called  the  zero  switch  rule  (ZSR)  (Maxwell,  24) 
and  has  generally  been  regarded  as  a  good  scheduling  policy.  The  ZSR 
seems  intuitively  to  be  sound  if  we  are  trying  to  minimize  inventory 
costs.  However,  if  we  find  an  optimal  solution  restricted  to  the  ZSR 
policy,  there  are  cases  where  we  can  improve  the  solution  by 
incorporating  a  non-zero  switch,  that  is,  for  some  product,  start 
production  before  its  inventory  reaches  zero. 

The  remainder  of  this  dissertation  examines  the  aforementioned 
issues  of  idle  time,  the  ZSR,  and  stochastic  input  to  a  bottleneck 
machine.  In  chapter  2,  a  review  of  relevent  literature  is  provided. 

This  includes  a  discussion  of  the  traditional  ELSP  and  a  review  of 
hierarchical  models  which  have  lot  scheduling  embedded  in  their 


structure. 
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Chapter  3  looks  at  two  important  issues  which  have  surfaced  in  the 
traditional  ELSP  literature,  idle  time  and  the  zero  switch  rule  (ZSR). 
In  particular,  we  find  that  the  conditions  where  inserting  idle  time 
improves  a  solution  to  the  ELSP,  given  zero  setup  costs,  are  very 
restrictive.  In  addition,  we  verify  that  the  ZSR  is  a  good  scheduling 
policy  for  most  problems  and  give  explicit  conditions  to  show  instances 
when  we  can  do  better  with  a  non-zero  switch. 

In  chapter  4,  we  examine  scheduling  a  bottleneck  machine  with 
stochastic  inputs  under  a  variety  of  situations.  First,  we  isolate  the 
issue  of  scheduling  deliveries  to  a  bottleneck  machine.  We  look  at 
using  current  state  information  to  schedule  the  pre-bottleneck 
machines.  Next,  we  consider  an  aggregate  planning  model  to  evaluate 
both  lot  scheduling  and  delivery  to  the  bottleneck  machine.  Then  we 
examine  a  dynamic  programming  formulation  under  varying  demand  and 
relax  the  constraint  that  all  demand  be  satisfied  during  each  period. 
Finally,  we  consider  a  variation  of  the  assembly  model. 

Chapter  5  develops  an  algorithm  for  solving  the  aggregate  planning 
model  of  chapter  4  and  shows  when  the  delivery  issue  is  important  as 
well  as  the  impact  on  the  system  of  variability  in  the  delivery 
process. 

Chapter  6  presents  conclusions  regarding  scheduling  idle  time,  the 
ZSR,  and  scheduling  a  bottleneck  machine  with  stochastic  inputs. 


CHAPTER  2 


LITERATURE  REVIEW 

2.1  The  Traditional  Economic  Lot  Scheduling  Problem  (ELSP) 

2.1.1  Background 

The  traditional  ELSP  addresses  scheduling  production  of  several 
products  on  a  single  machine.  Elmaghraby  (11)  provides  an  excellent 
review  of  the  traditional  ELSP  literature  through  1977.  The  following 
are  common  notation  and  assumptions  (with  i  being  the  index  referring 
to  a  particular  product): 

r^  demand  rate  in  parts  per  unit  time  (constant,  continuous) 

Pi  production  rate  in  parts  per  unit  time  (constant) 

Pi  =  ri  /  pi  relative  utilization  of  the  machine  by  product  i 

(  I  p.  $  1  for  feasibility) 

s^  setup  time  per  production  lot  of  product  i 

(assumed  independent  of  sequence) 

Aj  setup  cost  per  production  lot 

(assumed  independent  of  sequence) 
h^  holding  cost  per  part  per  unit  time 

T^  length  of  repeatable  production  cycle  for  product  i 

=  1/2  h^r^(l  -  p.)  scaled  holding  cost  when  T^  =  1 


The  average  cost  for  product  i  is 


Ci  =  Aj  /  Tj  +  hjrjd  -  PjjTj  /  2 

'  Aj  /  Tj  +  Hj  Tj  . 

The  objective  of  the  ELSP  is  to  find  T^  and  start  tines  for  each 
product  which  give  a  feasible  schedule  at  minimum  possible  cost.  A 
feasible  schedule  is  one  which  can  be  defined  on  a  Gantt  chart  over  any 
given  time  horizon  such  that  the  demand  for  each  part  is  satisfied 
throughout  the  time  horizon  and  the  machine  is  never  scheduled  for  more 
than  one  activity  at  any  point  in  the  time  horizon.  A  schedule  is 
infeasible  if  the  machine  is  required  to  work  on  more  than  one  product 
at  any  time. 

2.1.2  Independent  Solution  (IS)  to  the  ELSP 

Suppose  each  product  can  be  produced  on  independent  machines  with 
the  constraint  that  total  machine  time  used  is  equivalent  to  one 
machine.  Assume  inventory  of  each  product  follows  the  pattern 


One  Product  Inventory  Pattern 
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Average  inventory  cost  =  h^r^T.  (l-pp/2  =  ILTV  . 
Ve  can  formulate  this  problem  as 

(2-1)  minimize  £  H.TV  +  £  A./Ti 

(2-2)  subject  to  £  si/Ti  +  £  p.  £  1  . 


First  assume  that  (2-2)  is  not  binding  and  solve  (2-1)  to  get 

Ti*  =  (A./Hj)172  . 

If  this  solution  is  not  binding  in  constraint  (2-2),  then  the 
theoretical  optimum  solution  does  not  fully  utilize  machine  capacity. 
The  machine  capacity  utilized  in  this  solution  is 


0  =  X  SjOIj/Aj)1/2  +  X  Pi  . 


Observe  that  the  higher  we  set  the  values  of  A^,  the  less  machine 
capacity  we  use  in  the  theoretical  optimum  solution.  The  theoretical 
lower  bound  in  this  case  is 


C(Tt)  =  2  l  (A.IL)172  . 


Zero  setup  costs. 

In  the  case  of  zero  setup  costs,  use  Lagrangian  relaxation  on  (2-1)  and 


(2-2)  to  get 


T.  =  (s./H.) 

i  v  v  v 


{ l  (Vi>1/2} 

I1  -  i  >>} 
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Which  gives  the  lower  bound  on  solutions  of  this  form  with  zero  setup 
costs 

{ l  <Vi >1/2}2 

LLB  =  -  . 

H'*} 

2.1.3  Basic  Period  (BP)  Approach 

The  BP  approach  restricts  each  to  be  an  integer  multiple  of 
some  basic  period  w,  that  is,  =  K  w,  where  K  is  an  integer.  Several 
authors  have  developed  iterative  algorithms  using  this  approach,  for 
example,  Elmaghraby  (11),  Doll  and  Whybark  (10),  Fijita  (12),  and 
Madigan  (23).  These  iterative  techniques  involve  computing  a  basic 
period  w  and  rounding  cycle  times  T.  to  integer  multiples  of  w,  then 
checking  for  feasibility.  If  feasibility  is  not  achieved,  they  iterate 
to  a  new  basic  period  or  modify  the  solution.  A  BP  solution  is  feasible 
if  a  cyclic  production  schedule  can  be  found  such  that  the  requirements 
in  any  given  period,  including  setup  and  production  time,  do  not  exceed 
the  length  of  the  base  period  w.  Hsu  (11)  showed  that  the  check  for 
feasibility  is  NP-hard. 

2.1.4  Power-of-Two  Restriction 

Doll  and  Whybark  (10),  Delporte  and  Thomas  (7),  and  Fujita  (12) 
recommended  restricting  cycle  times  to  power-of-two  times  some  base 
period  to  simplify  the  search  for  a  feasible  schedule.  Maxwell  and 


Singh  (28)  showed  that  restricting  cycle  times  in  the  ELSP  to  power-of- 
two  times  some  base  period  would  yield  solutions  within  67.  of 
unrestricted  optimal  cycle  time  solutions.  Roundy  (29)  developed  an 
algorithm  to  perform  the  power-of-two  roundoff  of  order  intervals  which 
produced  solutions  within  6%  of  optimal  solutions.  Hence,  the  power-of- 
two  roundoff  technique  provides  a  simple  way  to  get  solutions  which  are 
very  close  to  optimality. 


2.2  ELSP  Imbedded  in  the  Hierarchy  of  Models 

The  key  feature  of  hierarchical  models  is  that  decisions  are  made 
at  one  level  under  consideration  of  their  impact  at  other  levels.  The 
linking  of  levels  of  decision  making  then  leads  to  better  overall 
solutions  to  the  system  being  modeled.  In  many  hierarchical  models 
involving  mathematical  programming,  this  linking  takes  the  form  of 
Lagrangian  multipliers  or  changes  to  the  constraints  (see,  for  example. 
Graves,  15).  Maxwell,  Muckstadt,  Thomas,  and  VanderEecken  (27)  propose 
a  general  modeling  framework  for  production  control  consisting  of  the 
following  three  phases:  creating  a  master  production  plan,  planning  for 
uncertainty,  and  real-time  resource  allocation.  Although  they  don't 
present  any  detailed  models,  their  approach  shows  how  various  models 
and  algorithms  can  effectively  be  placed  in  a  hierarchical  decision¬ 
making  structure.  Hax  and  Meal  (18)  give  reasons  for  using  a 
hierarchical  approach  and  develop  a  model  where  decisions  at  the 
aggregate  level  provide  constraints  for  the  lower  levels.  In  this 


model,  they  link  capacity  planning  decisions  with  detailed  scheduling 
decisions.  Grave's  model  (15)  decomposes  a  large  scale  production 
planning  problem  into  two  subproblems  corresponding  to  the  Hax-Meal 
model  and  also  provides  feedback  between  detailed  scheduling  decisions 
and  capacity  planning  decisions.  Dempster,  et  al.,  (8),  point  out  that 
the  two  fundamental  reasons  for  using  a  hierarchical  approach  are  to 
reduce  the  complexity  of  the  solution  process  and  to  cope  with 
uncertainty.  Bitran,  Haas  and  Hax  (2)  present  a  hierarchical  production 
planning  model  where  aggregate  planning  is  done  first  followed  by 
sequential  levels  of  disaggregation  in  the  production  planning  process. 
Maxwell  and  Muckstadt  (26)  develop  a  model  that  coordinates  production 
decisions,  including  capacity  planning  and  detailed  scheduling,  with 
transportation  decisions.  Hence,  we  see  that  the  hierarchical  approach 
provides  a  link  among  various  levels  of  decisions  in  a  manufacturing 
system. 


CHAPTER  3 


FURTHER  CONSIDERATIONS  ON  THE  TRADITIONAL  ECONOMIC  LOT 
SCHEDULING  (ELSP)  PROBLEM 

3.1  Then  to  Induce  Idle  Time  in  the  ELSP  with  Zero  Setup  Costs 

3.1.1  Introduction 
In  the  ELSP  with  zero  setup  costs,  intuition  suggests  that  machine 

capacity  should  be  fully  utilized  to  minimize  inventory  holding  costs. 
Dobson  (9)  presented  a  counterexample  in  which  an  optimal  solution  has 
idle  time.  This  section  addresses  the  conditions  under  which  it  is 
desirable  to  induce  idle  time  into  solutions  with  guaranteed 
feasibility.  Analysis  shows  that  these  conditions  can  generally  be 
characterized  as  one  product  having  dominant  holding  costs  and  the 
other  products  having  low  machine  utilization.  We  first  develop  the  two 
product  case  and  then  examine  two  production  patterns  of  the  three 
product  case.  These  lead  to  a  generalization  for  multi-products. 

3.1.2  Rationale  for  usins  zero  setup  costs. 

Since  the  fundamental  idea  of  the  ELSP  is  scheduling  multiple 
products  through  a  single  machine,  it  seems  only  natural  that  the 
dominant  feature  of  the  problem  is  the  constrained  resource  of  machine 
capacity  (see,  for  example,  Karmarkar,  21).  Therefore,  the  most 
important  impact  of  setup  arises  through  the  value  of  lost  machine 
capacity.  Hence,  in  this  section  setup  costs  are  set  to  zero  and  the 


constraint  on  machine  capacity  is  handled  explicitly. 


An  Aside  on  the  Baker  and  Bombereer  Problems 


The  classical  Baker  and  Bomberger  problems  (1,3)  have  high  setup 
costs,  that  is,  the  setup  costs  are  higher  than  the  value  of  machine 
capacity  lost  in  setup  time.  Hence  the  theoretical  optimum  solutions  to 
these  problems  have  substantial  idle  time  (see  Appendix  2).  Because  of 
this  idle  time,  several  authors  (7,10,11,13,14,23)  have  made 
significant  improvements  in  solving  these  problems  with  heuristic 
algorithms.  Clearly,  the  more  idle  time  you  induce  by  charging  high 
setup  costs,  the  easier  it  is  to  construct  feasible  schedules  that 
approach  the  theoretical  lower  bound. 


3.1.4  Two  product  case. 

Consider  two  products  produced  on  the  same  machine  where  one 
product  (call  it  product  1)  has  a  higher  holding  cost.  Suppose  for  a 
length  of  time  tQ  of  the  cycle  T,  we  produce  product  1  exactly  to  meet 
demand.  That  is,  idle  time  is  induced  after  producing  each  unit  of 
product  1  so  that  no  inventory  cost  is  incurred  on  product  1  during  the 
interval  (0,tQ)  (see  figure  3-1).  Conceptually,  you  could  also  view 
this  as  slowing  down  the  machine  during  tg  so  that  product  1  precisely 
meets  demand. 
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Product  1  |- 

Inventory  0 


FIGURE  3-1 

Two  Product  Inventory  Pattern 


The  time  parameters  tg,  t^,  t,,,  and  T  are  related  as  follows: 

t,,  =  ^ 2 ^  ^production  time  for  product  Z\ 

tj  =  [p1/(l-p1)](p2T+s1+Sg)  ^production  time  for  inventory 

cycle  of  product  l£ 

tf+Sg+tg+Si  =  (p^T+s^+s^)  /  (l-Pji)  ^length  of  inventory 

cycle  of  product  l£ 

t0  .  T  -  (t,+.2+t|8+»1) 

=  T(l-p1-pg)/(l-p1)  -  (s1+s2)/(l-p1)  . 

Since  tg  must  be  non-negative,  the  cycle  length  T  is  constrained  by: 


(1) 


T  >  (s,+s„)/(l-p1-pC))  . 


This  leads  to  the  overall  cost  expression  we  want  to  minimize: 

2 

Average  inventory  cost  =  H^(tj+32+t2+Sj)  ^  +  ^2^  or 

(2)  C(T)  =  H1|p22T+2p2(s1+s2)+(s1+s2)2/T|/(l-p1)2  +  HgT  . 

To  find  the  minimum  of  (2)  with  respect  to  T, 

(3)  dC/dT  =  H^/Cl-p^2  -  H^s^s^Vd-p^V  +  H2  . 

2  2 

Observe  that  d  C/dT  >  0,  hence  we  can  set  equation  (3)  equal  to  0  and 
* 

solve  for  T  to  find  the  minimum  of  C(T). 

T*  =  (s1+s2)  /  {(1-VVh,  *  p/}1/2 

Then  constraint  (1)  is  not  binding  and  tQ  >  0  in  the  optimal  solution  of 
(2)  onl:  f 

(Si+Sg)  /  |(l-Pj)^Hg/Hi  +  P22f2  >  (Sj+Sg)  /  (l-p^pg)  . 

This  gives  the  condition  for  inducing  idle  time 

H2/H1  <  (l-p1-2p2)/(l-p1)  . 

For  example,  if  p^  =  0.6  and  pg  =  0.15,  induce  idle  time  only  if 
Hg/Hj  <  1/4.  Observe  that  since  the  right  hand  side  of  this  inequality 
is  less  than  or  equal  to  1,  a  necessary  condition  for  inducing  idle 
time  for  product  1  is 


H0  <  H.  . 


f 
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Observe  also  that  since  the  left  hand  side  of  this  inequality  is  always 
greater  than  or  equal  to  0,  another  necessary  condition  is 

P2  <  (l-P1)/2  . 

Hence,  the  general  condition  for  inducing  idle  time  in  the  two  product 
case  is  when  the  product  with  lower  holding  costs  also  has  low  machine 
utilization. 

Observe  that  the  inventory  costs  using  induced  idle  time  may  be 

substantially  lower  than  the  theoretical  lower  bound  calculated  under 

the  assumption  of  full  saw-tooth  production  inventory  patterns. 

Consider  the  two  product  case  with  pg  <  (l-pp/2,  and 

Hg/Hj  <  (l-p1-2p2)/(l-p1).  Let  H2  =  H/k  ,  H1  =  kH  ,  therefore, 

2 

Hg/Hf  =  1/k  where  k  >  1  ,  sJ[  =  sg  =  s  ,  p1  -  0.5  ,  and  pg  <  0.25.  For 
this  range  of  parameters,  the  lower  bound  from  the  independent  solution 
is 


LLB 


|(Hs/k)1/2  +  (kHs)1/2j 

(l-.5-p  ) 

2 


=  gHsfk+n2  . 

k(l-p2) 

For  the  model  with  induced  idle  time, 

T*  =  2s  /  |(l-.5)2/k2  +  p22}1/2 

=  4k s  /  |l+4k2p22J1/2 


and 


.*  »-•  tJ  I.l  i.1  L*  L|  L|  Li  iL«  ViVI  mh 


C(T*)  =  8Hs|(l+4k2p22)1/2  +  2kp2J  . 


Then  we  can  compare  C(T*)  to  LLB. 

^ -  iS?  {(i+4tV)i/2  * 

If  we  take  the  case  when  p2  is  small,  we  get 

Liu  [C(t")/LLB]  «  4k  /  (k+1)2  . 

P2-  0 

Hence  in  this  case,  when  k  >  1  (that  is,  Hg  <  H^,  then  the  idle  time 
solution  is  less  than  the  independent  solution  lower  bound.  If  we  take 
the  case  when  k  is  large  (that  is,  Hg  «  H^),  we  get 


Lim  [C(T  ) /LLB]  =  16p2(l-p2)  . 


k  ® 


Hence  in  this  case,  when  Pg  <  .07,  the  idle  time  solution  is  less  than 

the  independent  solution  lower  bound.  If  we  take  the  example  where 

* 

P2  =  .05  and  k  =  5,  we  get  C(T  )/LLB  =  .85,  that  is,  the  average  cost 
per  year  for  the  idle  time  solution  is  less  than  the  corresponding  cost 
of  the  independent  solution  lower  bound. 
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In  summary,  we  see  that  in  the  two  product  case,  for  a  restricted 
range  of  the  parameters,  we  can  get  better  solutions  by  inserting  idle 
time  into  the  schedule.  The  restriction  on  the  range  of  the  parameters 
can  be  expressed  as  one  dominant  product  having  higher  holding  costs 
and  the  other  product  having  low  machine  utilization. 

3.1.5  Three  product  case. 

We  have  seen  for  the  two  product  case  conditions  under  which  it  is 
desirable  to  induce  idle  time  into  the  schedule.  Now  consider  three 
products  produced  on  the  same  machine  where  one  product  (product  1)  has 
higher  holding  costs.  There  are  an  infinite  number  of  possible 
production  patterns  we  could  consider.  However,  two  examples  of 
production  configurations  will  illustrate  that  the  conditions  for 
inducing  idle  time  in  the  three  product  case  are  similar  to  the  two 
product  case. 

3. 1.5.1  Production  Configuration  1. 

For  the  first  example  of  the  three  product  case,  we  consider  a 
production  pattern  with  two  setups  per  cycle  for  the  dominant  product. 
That  is,  we  build  up  inventory  of  product  1  prior  to  the  production  of 
each  of  the  other  products  (see  Figure  3-2). 
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The  time  parameters  tg,  t^2>  t^g,  tg,  tg. 


and  T  are  related  as  follows: 


t2  =  p2T 
*3  =  P3T 

t12  =  pl(p2T+8l+82)  7  (l-pP 
t12+S2+t2+Sl  =  (‘P2T+81+82^  7  (1_pP 

*13  =  pl(p3T+8l+83)  7  (1_pl) 

*13+83+t3+8l  =  ^p3T+81+83^  7  (1_pp 

tg  =  T(l-P1-P2-Pg)/(1-P1)  -  (2S1+S2+Sg)/(1-P1)  . 

Since  tQ  must  be  non-negative,  the  cycle  length  T  is  constrained  by 
(*)  T  £  (2s1+s2+8g)  /  (l-pj-pg-pg)  . 


Average  inventory  cost. 


H1  f  2  2  (s  +s2)2+(si+s'})2' 

C('T')  =  p  yZ  {T^P2  +p3  )+2P2(si+s2^+2p3^8l+83^  +  T 

+  HgT  +  HgT 

To  find  the  minimum  of  C  with  respect  to  T, 


(5)  dC/dT  =  H1|p22+Pg2-(s1+s2)2/T2-(s1+s3)2/T2|  +  H,,  +  Hg  . 


-2 


Observe  that  a  C/dT*  £  0.  Hence  we  can  set  equation  (5)  equal  to  0  and 
* 

solve  for  T  to  find  the  minimum  of  C(T). 


,+sJ  +(s4+sj‘ 


|(h2+h3)(i-p1)-/h1  +  Pg  +p 


Condition  for  inducing  idle  time 


Then  constraint  (4)  is  not  binding  and  tg  >  0  in  the  optimal 
solution  only  if 


{<s1+82)2+(s1+s3>211/2  2s.+s„+s„ 


|(h2+h3ki-p1)-/h1  +P8-+P3y'- 

Observe  that  (s1+s2)  +  (si+s3^  £  |^si+s2^2  +  ^sl+s3^2}1//2  * 

Furthermore,  note  that 

{(Hg+EjXl-Pj)*/*,  +  p/+P32}1/2  <  1-Pj-Pj-Pj 

only  if  (Hg+H^/H^  <  1.  Hence  a  necessary  condition  for  inducin 


time  is 


H2+^3  <  ^1  ’ 


Since  (Hg+Hg)/^  2  0,  another  necessary  condition  is 

k2+P  2)1/2  ,  k.+o2  ♦  (..♦o*)1'* 


s.+s«J+(s.+s 


V 


-22- 


If  we  take  the  simple  case  where  Pg  =  Pg  and  Sg  =  Sg,  this  reduces  to 

Pg+Pg  <  (l-p1)/2  . 

Note  that  this  is  similar  to  the  two  product  case.  The  general 
condition  for  inducing  idle  time  is  when  the  products  with  lower 
holding  costs  also  have  low  machine  utilization. 

3. 1.5.2  Dobson's  Counterexample 

Dobson  presented  the  following  counterexample  to  show  that  under 
certain  conditons,  it  was  optimal  to  have  a  schedule  with  idle  time 
(4,  page  18): 

hl  =  h2  =  h3  =  °*  rl  =  r2  =  r3  =  P1  P2  =  P3  =  S1  =  S3  =  i’ 

Sg  =  M,  p1  =  1/4,  p2  =  p3  =  1/M,  Hj  =  3/8,  and  Hg  =  Hg  =  0. 

These  parameters  clearly  satisfy  the  necessary  condition  (7), 
that  is,  0  +  0  <  1.  Idle  time  should  be  induced  if  constraint  (6)  is 
satisfied.  This  holds  for  values  of  M  as  small  as  5.  Hence,  we  don't 
need  to  use  extreme  pathological  cases  to  show  the  advantages  of 
inserting  idle  time. 

3. 1.5.3  Production  Configuration  2 

For  the  next  example  of  the  three  product  case,  we  consider  a 
production  pattern  where  the  dominant  product  is  only  set  up  once  per 
cycle.  That  is,  we  build  up  enough  inventory  in  that  one  production  run 
to  carry  through  the  production  time  of  the  remaining  products  (see 
figure  3-3). 


Product  1  |- 
Inventory  0 


Inventory  0 


FIGURE  3-3 


Three  Product  Inventory  Pattern 


„4  1L.*  U  '»J*1 


The  time  parameters  tQ>  t^  tg,  tg,  and  T  are  related  as  follows: 


l2  =  P2T 

l3  =  ^3T 

ll  =  Pi[(P2+p3)T+(si+s2+S3^  /  (1^) 

tl+s2+t2+s3+t3+sl  =  ((P2+P3)T+(S1+S2+S3)^  ^  (1-pi) 

*0  _  "^(l-Pi- Pg- P^)  ”  Cs^+Sg+Sg)/ Cl- • 

Since  tQ  must  be  non-negative,  the  cycle  length  T  is  constrained  by 


T  >  (Sj+Sg+Sg)  /  (l-prpg-p3) 


Average  inventory  cost. 


C(T)  *  Hi{T(P2+P3)2+2(p2+p3)+(8i'f82+S3)2/,T}  1  d'Pi)2  +  h2t  +  H3T 

To  find  the  minimum  of  C  with  respect  to  T, 

dC/dT  =  H1{(p2+p3)2-(s1+S2+s3)2/T2}  /  (1-Pl)2  +  Hg  +  H3  . 

2  2  * 

Observe  that  d  C/d  T  >0,  hence  we  can  solve  for  T  to  find  the 


minimum  of  C(T) 


(Si+S2+s3) 


{(H2+H3)(l-e1)2/H1  +  (P2+P3)2} 


2*1 1/2 


Constraint  (8)  is  not  binding  and  tg  >  0  in  the  optimal  solution 
only  if 

(H2+H3)/H1  <  1  -  2(p2+p3)/(1-p1)  . 

Observe  that  since  the  right  hand  side  of  this  constraint  is  less  than 
or  equal  to  1,  a  necessary  condition  for  inducing  idle  time  is 

H2+H3  <  H1  . 

Since  (H2+H3)/H^  >  0,  another  necessary  condition  is 

p2+P3  <  (l-pp/2  . 

Note  that  this  is  the  same  as  in  the  two  product  case.  Hence,  we  see 
that  for  general  production  configurations  of  the  three  product  case, 
the  two  necessary  conditions  for  inserting  idle  time  are  that  one 
product  has  dominant  holding  costs  and  that  the  remaining  products  have 
low  machine  utilization. 

3.1.6  General  criteria  for  inducing  idle  time. 

Using  the  results  derived  in  the  two  product  case  and  the  two 
production  configurations  considered  in  the  three  product  case,  we  can 
formulate  the  following  approach  to  determine  if  idle  time  should  be 
inserted  into  the  schedule  for  a  particular  problem: 


1.  One  product  (product  1)  has  dominant  holding  costs. 
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2.  Remaining  products  have  low  machine  utilization. 

I  Pi  <  (l-pp/2  • 

i  >  1 

If  a  problem  meets  these  necessary  conditions,  use  the  following 
procedure  to  determine  if  idle  time  should  be  induced: 

3.  Formulate  a  'good'  production  pattern  (for  example,  see  reference  6 
for  a  heuristic  technique  to  develop  production  sequences). 

4.  Define  the  relationships  between  time  parameters  t^  and  T  (similar  to 
those  defined  in  the  two  product  and  three  product  cases). 

5.  Formulate  the  average  inventory  cost  as  a  function  of  T  subject  to 

the  constraint  that  t^  is  non-negative. 

* 

6.  Solve  for  T  to  minimize  C(T). 

7.  If  the  parameters  of  the  problem  make  tQ  >  0,  then  idle  time  should 
be  induced. 

3.2  When  the  Zero  Switch  Rule  (ZSR-)  is  Not  Optimal 

3.2.1  Introduction. 

In  the  Economic  Lot  Scheduling  Problem  (ELSP) ,  the  ZSR  has 
generally  been  regarded  as  a  good  policy  for  keeping  average  inventory 
levels  low  (Maxwell,  24).  This  is  supported  by  the  fact  that  ZSR  is 
optimal  with  respect  to  inventory  of  a  single  product  as  will  be 
verified  in  the  following  section.  Pathological  cases  have  been 
developed  to  prove  the  non-optimality  of  the  ZSR  (Delporte,  6).  This 
section  addresses  specific  conditions  under  which  the  ZSR  is  non- 
optimal.  These  conditions  generally  occur  when  the  optimal  solution 
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with  respect  to  the  ZSR  gives  uneven  production  patterns  to  high  cost 
products. 


3.2.2  General  Approach. 

The  general  approach  we  use  to  improve  the  ZSR  solution  is  to  look 
at  marginal  adjustments  to  the  current  production  pattern  where  we  can 
decrease  the  overall  inventory  costs. 

3. 2. 2.1  Cost  Savings  from  Balancing  Production. 

Suppose  we  have  a  solution  to  the  ELSP  that  is  optimal  with 
respect  to  the  ZSR  policy.  Consider  two  adjacent  production  runs  of  one 
of  the  higher  cost  products  such  that  adjacent  runs  are  not  balanced. 


One  Product  Inventory  Pattern 


If  we  slightly  adjust  these  production  cycles,  retaining  the  same 
overall  length  T  for  both  cycles,  we  can  show  the  marginal  savings  in 
average  inventory  costs.  Consider  the  inventory  of  a  given  product 
under  the  ZSR  for  two  non-identical  production  runs. 


Inventory  f- - — 1 *f- - 

Level  T 

FIGURE  3-5 

One  Product  Inventory  Pattern 


Average  inventory  =  r(l-p) [t^2+(T-tj)2]  /  2T 

=  r(l-p)[2t12/T  +  T  -  2t^  /  2  . 

The  minimum  average  inventory  with  respect  to  t^  is  achieved  when  tj  = 
T/2,  that  is,  the  production  lots  are  balanced.  Then  the  cost  function 
can  be  written  as 

C  =  hr(l-p)[2t12/T  +  T  -  Zl  ]  /  2 

or  C(t1)  =  H[2t12/T  +  T  -  2t±]  . 

The  rate  of  change  of  C  with  respect  to  t^  is 

dC/dtj  =  4H[t1  -  T/2]  /  T  . 

Observe  that  if  t^  <  T/2  ,  increasing  t^  decreases  the  average 
inventory  and  decreasing  t^  increases  the  average  inventory. 

Conversely,  if  t^  >  T/2  ,  decreasing  t^  decreases  the  average  inventory 
and  increasing  t^  increases  the  average  inventory.  Hence,  balancing  the 


production  cycles  decreases  the  average  inventory. 


Consider  two  adjacent  production  runs  of  one  of  the  products  with 
lower  holding  costs.  Suppose  we  again  alter  the  cycle  to  accommodate  a 


non-zero  switch  while  retaining  the  same  overall  length  T  for  both 
cycles. 


Inventory  f 
Level  0 


t  t  1 

ll  1 


V'2 

Non-zero  switch 


FIGURE  3-6 


One  Product  Inventory  Pattern 


If  we  fix  T,  increasing  t^  decreases  the  length  of  the  second 
production  run.  We  can  show  the  marginal  increase  in  average  inventory 
costs.  Suppose  inventory  for  a  given  product  has  the  following  pattern: 


Inventory  f 
Level  0 


FIGURE  3-7 


One  Product  Inventory  Pattern 
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The  time  parameters  ^»tg,tg,t4,  ”  are  re^a^e<^  as  follows: 

tl+t3  =  or  tg  =  pT  -  t^,  and 

t2+t4  =  (l-p)T  or  t4  =  (l-p)T  -  t2  . 

Average  inventory  =  j  (p-r)t^2/2  +  [(p-r)t j+(p-r)tj-rtg]tg/2 

+  [(p-r)t1-rt2+rt4]tg/2  +  rt42  j  /  T 

=  Pt1t2/’r  "  rt2  +  r(l-p)T/2  . 

Observe  that  t^  £  pT.  Hence  the  coefficient  of  tg  is  always  less  than 
or  equal  to  0.  Because  the  inventory  at  the  non-zero  switch  point  must 
be  non-negative,  we  get  the  following  constraint: 

(p-r)tjl  -  rt2  £  0  . 

This  constrains  tg  to 

0  £  t2  £  t1(p-r)/r  . 

Hence,  average  inventory  of  this  product  would  be  minimized  with 
respect  to  tg  by  tg  =  t^Cp-O/r.  That  is,  average  inventory  of  that 
particular  product  is  minimized  by  following  the  zero  switch  rule.  Note 
also  that  average  inventory  of  this  product  is  then  minimized  with 
respect  to  tg  (under  the  ZSR)  when 

=  pT/2  . 


-3J 


That  is,  average  inventory  is  minimized  when  the  production  lots  are 
exactly  balanced.  Then  the  cost  function  can  be  written  as 

C  =  h^ptjtg/T  -  rtg  +  rT(l-p)/2*  . 

The  rate  of  change  of  C  with  respect  to  t^  is 

dC/dt1  =  hptg/T  . 

3. 2. 2.3  Improving  the  ZSR  Solution. 

We  can  improve  the  ZSR  solution  in  the  following  way.  Take  two 
products  that  have  two  adjacent  production  runs  in  the  ZSR  solution. 
Adjust  their  production  times  retaining  ZSR  for  the  high  cost  product 
(product  i)  and  incorporating  a  non-zero  switch  for  the  lower  cost 
product  (product  j),  as  shown  in  Figures  3-4  and  3-6,  such  that  we 
don't  affect  the  rest  of  the  production  cycle.  We  then  make  the 
production  pattern  of  the  higher  cost  product  more  balanced  at  the 
expense  of  using  a  non-zero  switch  policy  for  the  lower  cost  product. 

Then  we  should  use  the  non-zero  switch  option  if  the  marginal  cost 
savings  from  making  production  of  product  i  more  balanced  exceeds  the 
marginal  cost  of  using  a  non-zero  switch  on  product  j  .  That  is 


2H.(T.-2t..)/T.  >  h.p.t.0/T. 

iv  l  ir '  l  j*j  j2'  j 


or  (9) 


Hi  2ri(1-Pi>  ,  Ti  Hj2 

*4  Pj  ITj/S-ti  -p  Tj 


t: 


Thus  we  see  that  we  can  improve  a  ZSR  solution  when  that  solution  gives 
a  lumpy  production  pattern  to  a  product  with  dominant  holding  costs. 

By  using  a  non-zero  switch  on  some  cheaper  product,  we  can  balance 
production  of  the  dominant  product. 

3.2, 2. 4  Delnorte's  Counterexample. 

Delporte  presented  a  counterexample  to  the  optimality  of  the  ZSR 
(6,  App  IV,  A  Counterexample  to  the  Optimality  of  the  ZSR).  The  data 
for  this  problem  is  as  follows:  pg  »  1,  rg  ■  .2,  pg  =  .2,  pg  =  1, 
r3  =  **»  P3  =  •*»  1*3  =  ♦!»  ri  ~  =  *4,  and  h^  =  100.  The  solution 

to  this  problem  with  respect  to  the  ZSR  is  given  by  Figure  3-8. 
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0  4  6.462  10.155  15.693  20 

FIGURE  3-8 

Three  Product  Inventory  Pattern 


Substituting  this  data  into  equation  (9),  we  get  480  >  1.3.  In  this 
example,  we  see  that  the  ZSR  solution  is  not  optimal  even  if  we  take 
less  extreme  values  of  relative  holding  costs,  that  is,  as  long  as 
hi/h3  >  370  ,  the  ZSR  solution  is  not  optimal. 

3.2.3  Conditions  When  the  ZSR  is  not  Optimal 
Under  what  conditions  is  it  likely  that  the  ZSR  is  not  optimal?  We  can 
best  answer  that  by  asking  when  the  inequality  (9)  is  likely  to  hold. 
Consider  two  products  whose  independent  solutions  for  natural  cycle 
length  are  approximately  equal,  that  is,  using  the  Lagrangian 
relaxation  method,  we  obtain 

Ti'  ■  { l  t=i/»i31/s  /  {i-X  Pi}  • 

Hence,  these  two  products  should  be  produced  at  roughly  the  same 

*  * 

relative  frequency.  Then  T.  =  T.  or  s./H.  =  s./H..  If  we  let 

l  j  i  i  J  J 

Tj  =  Tj  and  substitute  into  equation  (9),  we  get 

(,./s.)  p.(l-p.)  >  tj2/(T/2  -  tu)  . 
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Then  under  this  situation,  the  ZSR  solution  is  not  optimal  if 

(a)  Two  products  have  roughly  the  same  natural  cycle  length 

while  one  of  the  products  has  both  dominant  holding  costs 
and  setup  time. 


(b)  The  ZSR  solution  gives  uneven  production  patterns  to  the 
higher  cost  product. 


In  summary,  conditions  (a)  and  (b)  above  describe  when  the  ZSR  is 
likely  to  not  be  optimal  and  we  may  be  able  to  improve  the  solution  by 
incorporating  a  non-zero  switch.  When  these  conditions  do  not  hold,  we 
can  safely  conclude  that  the  ZSR  is  a  good  policy.  For  example,  if  we 
use  a  power-of-two  policy  (Roundy,  29),  all  production  lots  for  a  given 
product  are  equal  so  we  could  not  improve  the  solution  by  using  a  non¬ 
zero  switch. 


CHAPTER  4 


STOCHASTIC  INPUT  TO  A  BOTTLENECK  MACHINE  i 

i 

4.1  Introduction  j 

! 

Consider  a  bottleneck  machine  which  processes  several  products  | 

(see  Figure  4-1).  The  traditional  ELSP  approach  is  to  assign  costs  for  i 

I 

setup  of  each  product  and  costs  for  holding  inventory  after  processing  j 

until  the  parts  are  demanded  or  consumed.  It  ignores  the  issue  of 

delivery  of  raw  parts  to  the  bottleneck  machine  and  resulting  holding 

costs  for  those  raw  parts.  If  delivery  of  the  raw  parts  is  j 

deterministic,  the  delivery  of  each  product  could  be  scheduled  to 

coincide  exactly  with  its  start  time.  In  this  deterministic  case,  it  is 

reasonable  to  ignore  the  delivery  issue  when  scheduling  the  bottleneck 

machine.  However,  if  the  deliveries  are  not  deterministic,  that  is  if 

there  is  variability  in  shipping  time  or  processing  time  on  the 

predecessor  machines,  then  ignoring  this  issue  when  scheduling  the 

! 

bottleneck  machine  can  seriously  affect  total  operating  costs.  For 
example,  if  delivery  of  the  raw  parts  is  requested  too  early,  then 
excessive  work-in-process  (IIP)  inventory  will  accumulate  before  the 
bottleneck  machine.  On  the  other  hand,  if  delivery  of  the  raw  parts  is 
requested  too  late,  then  productive  capacity  of  the  entire  system  is 
diminished  because  the  bottleneck  machine  is  delayed. 

This  chapter  addresses  the  combined  issues  of  scheduling  delivery 
of  raw  parts  to  a  bottleneck  machine  and  scheduling  lot  sizes  on  that 


machine.  Section  4.2  focuses  on  the  delivery  issue  assuming  the 
production  schedule  has  already  been  defined  for  the  bottleneck 
machine.  Section  4.3  analyzes  how  state  information  can  be  used  to 
schedule  the  machines  using  real  time  data.  Section  4.4  develops  an 
aggregate  model  for  combining  the  issues  of  lot  sizing  and  delivery 
scheduling.  Section  4.5  relaxes  the  constraint  that  all  demand  be 
satisfied  in  each  period  and  looks  at  a  dynamic  programming  approach 
Section  4.6  looks  at  the  special  case  of  an  assembly  model  with  a 
different  inventory  pattern. 

Predecessor  Machines 


FIGURE  4  -  1 


Machine  Network 
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4.2  Scheduling  Deliveries  to  a  Bottleneck  Machine 


4.2.1  Simple  Newsboy  Model 

Consider  one  particular  production  lot  on  the  bottleneck  machine 
and  focus  only  on  the  delivery  issue.  Suppose  the  machine  will  be 
available  and  production  of  this  lot  is  scheduled  to  start  at  a 
specific  time  (reference  this  as  time  zero).  Define  the  following 
additional  parameters: 


requested  delivery  time  of  product  i  raw  parts  , 

random  variable  representing  actual  delivery  time 
(note:  D^  >  0  means  the  lot  arrived  late,  <  0 
means  the  lot  arrived  early)  , 

probability  density  of  delivery  time  of  product  i 
centered  about  mean  zero  (assumed  independent  of  lot 
size  —  this  assumption  may  not  be  valid  in 
certain  cases,  Karmarkar,  22)  , 
cumulative  distribution  of  f^(.)  , 

holding  cost  at  the  bottleneck  machine  of  raw  part  i 

per  unit  time  , 

lot  size  for  product  i  ,  and 

value  of  lost  machine  capacity  per  unit  time  . 


V-> 


pi<-> 


Given  that  d^  is  the  requested  delivery  time,  there  are  two  possible 
outcomes  and  costs  incurred. 


> 

t 


Case  1.  Early  delivery  (D.^  <  0) 


The  cost  incurred  is  h (— D^)  which  is  the  inventory  holding  cost 
of  raw  parts. 

Case  2.  Late  delivery  (D^  >  0) 

- T - T - > 

0  t 


The  cost  incurred  is  X  which  is  the  value  of  lost  machine  capacity. 
The  expected  cost,  given  that  d^  is  the  requested  delivery  time,  can  be 
formulated  as 


ru  r® 

C(d.)  =  h.r.T.  -x  f . (x  -  d.)  dx  +  X  x  f.(x  -  d.)  dx  . 

^  ^  *  X  x  X  ,  ft  X  1 

—CD  m  0 


Using  a  change  of  variable,  y  =  x  -  d^  ,  we  simplify  this  expression 
to 


d .  CD 

C(di)  =  -  h.r.r  j  1  (y  +  d.)  f ,(y)  dy  +  X  J  (y  +  d.)  f .(y) 

i 


dy 


_  r  dj  p  d. 

=  X  d.  -  (X  +  h.r.T.)  y  f.(y)  dy  -  (X  +  h.r.T.)  d.  |  f  (y)  dy  . 

0  —CO  J  —co 


To  find  the  optimal  requested  delivery  time,  determine 


Since 


d2  C 


ddi 


o  =  (X  +  h.r.T.)  f.(-d.)  £  0  ,  C  is  convex  in  d.. 

m  11111  1 


Then  the  optimal  requested  delivery  time  is  to  choose  d^  such  that 


Fi(-di  >  - 


X  +  h.r.T. 

Ill 


This  is  the  familiar  solution  to  the  newsboy  problem.  In  this  case,  we 
are  balancing  the  expected  value  of  lost  machine  capacity  with  the 
expected  value  of  inventory  holding  costs. 


4.2.2  Stochastic  Delivery  with  Two  Random  Variables 

Consider  one  production  lot  (product  i)  scheduled  to  start  at  a 
specific  time,  reference  this  as  time  zero,  and  another  production  lot 
(product  j)  scheduled  to  start  upon  completion  of  product  i  (see 
figure  4-2). 


Scheduled 


Actual 


Figure  4-2 
Gantt  Chart 


If  <,  0,  the  production  of  lot  i  will  start  on  time,  that  is,  time 
zero,  and  if  >  0,  it  will  start  late.  Conditioning  on  D^,  there  are 
again  two  possible  outcomes  and  costs  incurred  due  to  product  j. 


Case  1.  Early  Delivery,  that  is,  <  Pi  +  max?0,  D.f  . 

The  cost  incurred  is  -  h.r.T.(D.  -  P.  -  max$0,  D.f)  . 

JJJ  J  1  * 

Case  2.  Late  Delivery,  that  is,  D^  >  +  max$0,  D^|  . 

The  cost  incurred  is  X  (D^  -  Pi  -  max$0,  D^O  . 


The  expected  cost,  given  that  d^  is  the  requested  delivery  time  and  f\ 
is  the  distribution  of  the  start  time  of  production  lot  i,  can  be 
calculated  as 


vv  ■  J0  riw  dy  [  -v^  fS y  <x ' Pr  y)  V 


V  ** 


CO 

f  (x  -  P  -  y)  f  (x  -  d.)  dx  ]  . 

Jp  +  y  1  J  J 


Using  a  change  of  variable,  z  =  x  -  d.,  this  simplifies  to 

J 


W  *  J0  fi(y)  dy  [  -KjrjTj 


y  -  d 

J  (z  +  d^  Pi~  y)  f j  (z)  dz 


'.+  y  -  d . 
i  J 


(z  +  d  -  P  -  y)  f  (z)  dz  ]  . 


The  optimal  requested  delivery  time  d.  is  found  using  the  following 

J 

equation: 


r 

=  X  -  (X  +  h  .r  .T .)  f .  (y  )  F,(P;+  y  -  dj  dy  =  0 


d  C. 
J 


[>■< 
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I f  f\  represents  the  arbitrary  start  time  of  production  lot  i, 
incorporating  all  random  events  up  to  that  point,  then  the  same  form  of 
solution  applies. 


4.2.3  Three  Echelon  Machine  Network 

Suppose  the  inputs  to  the  predecessors  of  the  bottleneck  are 
stochastic  (see  figure  4-3).  The  issue  is  when  to  schedule  deliveries 
to  machine  i.  If  the  delivery  arrives  early,  IIP  is  added  to  the  system 
at  machine  i.  If  the  delivery  arrives  late,  the  subsequent  delivery 
from  i  to  the  bottleneck  might  be  late  thereby  causing  a  reduction  in 
capacity  of  the  overall  system.  Since  we  know  the  rate  of  change  of 
cost  at  level  i  with  respect  to  d^,  this  is  the  marginal  rate  at  which 
the  overall  system  costs  would  be  affected  by  a  delay  in  delivery  to  i. 


Bottleneck 

Machine 


FIGURE  4-3 


Machine  Network 


The  overall  expected  system  cost  with  respect  to  delivery  from  machine 
ij  can  be  expressed  as 


cu<du>  -  J_ 1  Vx  -  V  dI 


00 


[X  -  (X  +  h.r.T.)F.(-d.)]  f  x  f .  .(x  -  d.  .) 

iiii  i  J0  ijv  ij' 


.  dx 


=  [X  -  (X  +  h.r .T.)  F.(-d.)  ]  d.  . 
L  i  i  i  iv  r  J  ij 


„-d. 


-  +  Kijri jTi j  -  <X  +  ViV  V'V  1  JJj  *  fijW 


dy 


N-d .  . 


-  +  VijTij +  (x  +  EiriTi)  Fi(-di>  i  JJj  ^  ■ 


This  yields  the  optimal  delivery  time  d^  that  satisfies 
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In  summary,  delivery  times  should  be  chosen  to  balance  the  expected 
costs  of  holding  inventory  with  the  expected  cost  of  lost  machine 
capacity.  Having  completed  this  analysis  for  computing  deliveries  for 
general  points  in  time,  we  now  develop  a  method  that  looks  at  the 
current  state  of  the  machine  network  for  real  time  scheduling  of 
deliveries. 

4.3  Using  State  Information  to  Schedule  Pre-Bottleneck  Machines 

4.3.1  General  Formulation 

Using  the  concepts  developed  in  Section  4.2,  we  can  then  use 
current  information  on  the  state  of  the  system  to  schedule  production 
in  the  predecessors  to  the  bottleneck  machine.  Suppose  a  given  machine 
i  has  finished  production  of  a  lot  and  sent  it  to  the  bottleneck 
machine.  When  should  we  initiate  a  new  production  lot  of  i  knowing 
current  state  information,  i.e.  the  status  of  the  bottleneck  machine 
and  other  predecessors  to  the  bottleneck?  In  general,  the  expected 
delivery  time  of  the  raw  parts  should  match  the  expected  availability 
of  the  bottleneck  machine  plus  an  allowance  for  some  safety  time.  The 
amount  of  safety  time  depends  on  the  relative  value  of  machine  capacity 
to  the  cost  of  holding  parts. 

Let  g^(x  |  state)  be  the  conditional  distribution  of  the  time  the 
bottleneck  machine  will  be  available  again  for  product  i,  given  the 
current  state  of  the  system.  Then  the  contribution  of  the  delivery  of 
product  i  to  the  system  cost  can  be  formulated  as 


Ci(di)  =  J0  gi^x  I  8tate)  C-^iriTi  J  (y-*)  fi(y-di)  dy 


r® 

+  x  (y-x)  f^y-d.)  dy  ]  dx  . 
J  y 


The  optimal  solution  is  to  choose  cL  such  that 


03 

j  g.(x  |  state)  F . (x-d.)  dx  =  - - - 

JO  1  11  X  +  h.r.T. 


ill 


To  shorten  future  notation,  let  R^  = 


X  +  h.r.T. 

ill 


4.3.2  Two  Product  Case 

A  two  product  case  illustrates  how  current  state  information  can 
be  used  to  schedule  pre-bottleneck  machines.  Consider  two  machines  A 
and  B  which  directly  precede  the  bottleneck  (see  figure  4-4). 


FIGURE  4-4 
Machine  Network 


Suppose  machine  A  has  finished  and  sent  a  production  lot  to  the 
bottleneck  (see  figure  4-5).  When  do  we  initiate  production  of  A 


< - > 

< - 1 - > 

<1 -  PA - >< - PB - > 

decision  point  available  for  A 

FIGURE  4-5 
Gantt  Chart 

Ve  want  to  know  how  much  work  remains  on  the  bottleneck,  including  the 
distribution  of  possible  delays  as  a  function  of  state  variables, 
before  the  bottleneck  will  be  available  for  product  A  again. 

The  possible  states  for  A  and  B  can  be  defined  as  follows: 

Production  lot  A. 

State  1.  Finished  on  the  bottleneck. 

State  2.  Started  but  not  finished  on  the  bottleneck. 

State  3.  Not  started  on  the  bottleneck. 

Production  lot  B. 

State  1.  Started  on  the  bottleneck. 

State  2.  Awaiting  processing  on  the  bottleneck. 

State  3.  Started  but  not  finished  on  machine  B. 

State  4.  Not  started  on  machine  B. 

The  state  of  the  system  can  be  expressed  as  (a,b)  where  a  represents 
the  status  of  production  lot  A,  and  b  represents  the  status  of 
production  lot  B.  Observe  that  we  don't  need  to  consider  state  (1,2) 
because  B  has  nothing  to  wait  for  in  this  case  and  states  (2,1)  and 
(3,1)  because  of  incompatibility.  For  each  of  the  possible  states,  we 
can  then  formulate  an  expression  corresponding  to  the  optimal  requested 
delivery  time  developed  in  Section  4.4.1.  For  each  of  these  states,  we 


Machine  A 
Machine  B 

Bottleneck 
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will  illustrate  the  results  with  exponential  and  uniform  distributions 
of  processing  times  on  machines  A  and  B  with  constant  deterministic 
processing  times  on  the  bottleneck  machine. 


State  (1.1).  A  has  finished  on  the  bottleneck,  and  B  has  started  on 
the  bottleneck. 


Let  t^  be  the  time  B  has  already  been  processed  on  the  bottleneck 


Bottleneck  <- 


T 

t. 


available  for  A 


The  optimal  solution  is  to  choose  d  such  that  F  (P„-  t.  -  d  )  =  R 

a  a  B  b  a  a 


Case  1.  Exponential  processing  times  for  machines  A  and  B. 

Assume  A  and  B  have  exponential  processing  times  on  the  pre-bottleneck 
machines,  with  deterministic  processing  times  on  the  bottleneck 
machine,  that  is,  for  machines  A  and  B, 


f^(x)  =  e  Xix  ,  for  x  £  0  , 


and  F^x)  =  1  -  e  Xix  ,  for  x  £  0  ,  i  =  A  and  B. 


The  optimal  solution  in  this  case  is  to  choose  d  such  that 

a 


1  -  e~Xft(PB  dft)  =  R 


or 


d  *  -  P»  “  t.  +  i-  log  (1  -  R.)  . 


Since  0  £  R  £  1  ,  log(l  -  R  )  £  0  .  Observe  that  this  corresponds  to 

&  a 

matching  the  delivery  time  with  the  availability  of  the  bottleneck 

machine  (in  this  case,  Pg  -  t^)  plus  allowing  for  allowing  for  safety 

time  (in  this  case,  l/X  log  (l  -  R  )). 

&  & 


Case  2.  Uniform  processing  times  for  machines  A  and  B 


For  machines  A  and  B, 

f^x)  =  1  /  2u.  ,  for  0  <,  x  £  2ua  , 

and  F^(x)  =  x  /  2u^  .for  0  £  x  £  2u^  . 

The  optimal  solution  in  this  case  is  to  choose  d  such  that 

d 

<V  V  d.)  /  2\  - 

* 

or  d  =  PD  -  t,  -  2u  R  . 
a  B  b  a  a 

Observe  that  this  corresponds  to  matching  the  delivery  time  with  the 

availability  of  the  bottleneck  machine  (in  this  case,  Pg-  tb>  plus 

allowing  for  safety  time  (in  this  case,  2u  R  ). 

a  a 


State  (1,3).  A  has  finished  on  the  bottleneck,  and  B  has  started  on 


machine  B. 

Let  tg  be  the  time  B  has  already  been  processed  on  machine  B  and 

fg(x  (  tg)  be  the  distribution  of  time  remaining  on  machine  B. 
Machine  B  < - 1 - > 


Bottleneck 


available  for  A 
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The  optimal  solution  is  to  choose  d^  such  that 

J0  V*  I  tB)  P,<x  +  V  V  *  R*  • 

Case  1.  Exponential  processing  times  for  machines  A  and  B. 
The  optimal  solution  is  to  choose  d&  such  that 


[i  - 


-Xa<x  +  pB- 


t,  -  d  ) 
b  a7 


]  dx 


R 


a 


°r  d»‘  -  pb  +  -r  +  -f  l°« -  p»>  • 

aba 


Case  2.  Uniform  processing  times  for  machines  A  and  B. 

Assume  that  the  distribution  fg(x  |  t0)  is  uniform  (0,  2ub~2tg) 

The  optimal  solution  is  to  choose  dft  such  that 


f 

Jo 


2ub"2tB  1  rX  *  PB~  da1  =  p 

2u, -2tp  2u  J  a 

b  a  a 


or  da  *  PB  +  V  lB  -  2u«B.  ’ 


State  (1.4).  A  has  finished  on  the  bottleneck,  and  B  has  not  started 
on  machine  B. 

Let  fgjC*)  be  the  distribution  of  start  time  on  machine  B. 

Machine  B  < —  x  — >< —  y  - > 

Bottleneck  < -  P,, - > 


available  for  A 


The  optimal  solution  is  to  choose  d&  such  that 


-OB  .« 

Jo  f»1<x>  Jo 


fB(y)  pa(x  +  y  +  Pp-  d.)  dy  dx  =  R.  • 


B  ua' 


i 
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Case  1.  Exponential  processing  times  for  machines  A  and  B. 


The  optimal  solution  is  to  choose  d&  such  that 


_xbix  r*  .  _xsy 


i 

e  __  X^e  D  [1  -  e 

Jo  15 


-y* +  y +  V  V  V 


]  dy  dx  =  R, 


or  d»  =  Pb  +  “T"  log[-  lost 


X-+Xbi  ,  1  ,._rXa+Xbl-,  ,  1 


^3 


y  ]  +  ~ '«  (1  -  v- 


Case  2.  Uniform  processing  times  for  machines  A  and  B. 

Assume  that  the  distribution  fg(x  |  tg)  is  uniform  (0,  2ub~2tg) . 


The  optimal  solution  is  to  choose  d  such  that 

a 


r2ubi  _j_  r2ub  j_  r 

Jq  2uhi  Jn  2uh 


.2u, 

bl  Jo 


1  +  1  4  PB-  d. 


Oil 


]  dy  dx  =  R 


or  d>  =  PB  +  "b  +  “bl  '  2“»»,  ' 


State  (2.2),  A  has  started  on  the  bottleneck,  and  B  is  awaiting 
processing  on  the  bottleneck. 


Bottleneck  <- 


t 


-><- 


B 


available  for  A 


The  optimal  solution  is  to  choose  d  such  that 

a 


F  CP.-  t  +  PB-  d  )  =  R  . 

a'  A  a  B  a'  a 


Case  1.  Exponential  processing  times  on  machines  A  and  B. 


The  optimal  solution  is  to  choose  d  such  that 

a 


1  _  e-y  v  v  v  y  .  B 
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or 


d,'  *  PA  -  ‘a  +  PB  +  4-  l0*  -  V 

a 


Case  2.  Uniform  processing  times  on  machines  A  and  B. 


The  optimal  solution  is  to  choose  d  such  that 

Si 


(P.-  t  +  PD-  d  )  /  2u  =  R 
v  A  a  B  a'  '  a  a 


or 


d  =  P.  -  t  +  PD-  2u  R  . 
a  A  a  B  a  a 


State  (2. O’).  A  has  started  on  the  bottleneck,  and  B  has  started  on 
machine  B. 

Machine  B  < - ■ - > 


Bottleneck 


t" 


-X- 


B 


-> 


a 


available  for  A 


Then  the  distribution  of  the  start  time  of  B  on  the  bottleneck  is 
defined  by  the  maximum  of  the  completion  time  of  B  on  machine  B  and  the 
completion  time  of  A  on  the  bottleneck. 


Let  fgb(x  I  ta,th)  be  the  distribution  of  the  start  time  of  B  on  the 


'a’  b 


bottleneck  and  fgb(x  |  ta>tb)  =  <  Fb(x  |  tb)  for  x  =  P^-  t 

fb^x  I  for  x  >  PA~  la 

The  optimal  solution  is  to  choose  d  such  that 


W  *.  I  V  W  V  pB-  V 

+  Jp  _t  fb(X  1  V  Fa(*  +  PB-  d.)  dx  - 


-GO 


A  a 


-5 


Case  1.  Exponential  processing  times  on  machines  A  and  B. 
In  the  exponential  case,  we  get 


I  V  v  - 


1  -  e 


,  for  x  =  P . -  t 
A  a 


,  for  x  >  P.-  t 
A  a 


The  optimal  solution  is  to  choose  d&  such  that 

r  “  VPA_  V  n  _Xa(PA~  V  PB  "  da) 

[1  -  e  J  [1  -  e  J 


L  t  v 


[1  -  e 


-Xa(x  +  Pg  -  d&) 


]  dx  =  R 

J  • a 


A  a 


d  =  P.-  t  +  Pn  -  -i-  log[l  -  r-4j-  e 
a  A  a  B  X&  VS 


Xa  "VPA  V-,  1  .  n  p  v 

e  ]  +  “^lofi  (1-Ra) 


Case  2.  Uniform  processing  times  on  machines  A  and  B. 

In  the  uniform  case,  we  get 

fsb^X  I  W  =  2u,  -  2t,  ’  f0r  X  =  PA~  ta  * 

D  D 

2u,  -  2t.  ’  for  Pa"  ta  <  x  *  2ub~2tb  • 

b  b 

The  optimal  solution  is  to  choose  d  such  that 

a 


V  la  rVV  PB~  da 
2V2tBL  2ua 


r2u,  - 

]  +  J 

JP,-t 


,2u, -2t.  ,  x  +  Pn-  d 

P.-t  2V2tbC  2ua  ]  ‘  ! 

A  a 


°r  d.  *  PB  +  Ub'  VTOTCT-2“.Ra 


State  (2.4).  A  has  started  on  the  bottleneck,  and  B  has  not  started 


°r  ■*/  -  pA  +  pb  +  r  l0‘  -  V  +  r  l0‘  <^7^ 

a  a  al 


Case  2.  Uniform  processing  times  on  machines  A  and  B. 


The  optimal  solution  is  to  choose  d&  such  that 


r2ual  1 
•*0  2ual 


(x  +  PA+  PB-  da) 


or  d  =  P.  +  P_  +  u  «-  2u  R 
a  A  B  al  a  a 


State  (3,3).  A  has  not  started  on  the  bottleneck,  and  B  has  started 


machine  B. 


Machine  B  <- 


Bottleneck 


< —  x  - ><- 


PA - ><- 


available  for  A 


fgb(k)  =  the  distribution  of  max$x  +  PA,  y\. 


Case  1.  Exponential  processing  times  on  machines  A  and  B. 


In  the  exponential  case,  we  get  for  k  ^  P, 


V  rk  *  Xal(x  V 


*  V  ^  I  ^ 

"a 


Xai^k  pa^  rk .  “V 


,  .  A  alv  ^  r .  'V  . 

+  Xale  J  V  dx 


V  .  ,  .  Xal(k  “  PA)  (  x  _Xbk  "  Xal(k 

-  (Xal+Xb)  e 


=  V  +  Xale 


J2 


■  V  ’/  V  V  .* 


The  optimal  solution  is  to  choose  d  such  that 

a 


f  fsboo  [1  -  . 


-xa(k  +  pb-  da) 


]  dk  =  R. 


or 


da  =P.+P 


X  XX  _XbPA 

*+  P„-  log[r — - 77  ,a*  y7T~TT "  j.'  \  ~\]  +  7“  log(l-R  ). 

A  B  Xa  Xal+  Xa  ^VV^VV  Xa 


Case  2.  Uniform  processing  times  on  machines  A  and  B. 

The  optimal  solution  is  to  choose  d  such  that 

a 


I 


,2u.  -  t,  . 
b  b  r  k 


L2u  ,  2u. -2t,  -  P. 
al  b  b  A 


2u,  -  t,  2u  < 
b  b  al 


,  k  +  PR-  da 

■]  c— W  ■■  ] dk 


f,PA+  2ual  _ 1_  rk  *  PI 

•*2u.  -  2t,  2ual  2us 


V  da. 


=  R 


or  d  *  =  P  -  4UalUaRa  +  3(2V  tb)PA 

8  B  2\r  pa2/  <*v*V  V 


State  (3,4).  A  has  not  started  on  the  bottleneck,  and  B  has  not 
started  on  machine  B. 

Machine  B  < —  z  — >< -  y  - > 

Bottleneck  < - x  - >< -  P^ - >< -  Pg  - > 

available  for  A 

fgb(k)  =  the  distribution  of  max$x  +  P ,  y  +  z\ 


The  optimal  solution  is  to  choose  d  such  that 
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Case  1.  Exponential  processing  times  on  machines  A  and  B. 

Similar  to  state  (3,3)  except  X^  is  replaced  by  ^  +  X^. 

Case  2.  Uniform  processing  times  on  machines  A  and  B. 

Similar  to  state  (3,3)  except  with  added  time  delay  at  machine  B. 
In  general,  the  use  of  state  information  to  determine  delivery  time  can 
be  characterized  as  matching  the  expected  delivery  time  to  the 
availability  of  the  bottleneck  machine  in  addition  to  allowing  for 
safety  time  dependent  on  the  relative  value  of  machine  capacity  to 
holding  cost. 

4.4  Aggregate  Planning  Model  with  Lot  Sizes  and  Deliveries  as 
Variables. 

In  the  previous  sections,  delivery  times  were  determined  when  the 
lot  sizes  were  known.  In  this  section,  both  delivery  times  and  lot 
sizes  will  be  considered  as  variables.  The  traditional  ELSP  (relaxing 
scheduling  constraints)  is  formulated  as  follows: 

minimize  7  (H.T.  +  A./  T.) 

L  v  i  i  r  v 

subject  to  ( s i /  Tj  +  p.)  £  1  . 

In  this  model,  the  optimal  lot  size  is  given  by  rp^.  The  average 
inventory  cost  when  delivery  times  are  random  variables  is 

1  -  r° 

l  —  h.rp.  -x  f.(x  -  dp  dx  . 

i  “  —o> 

To  shorten  future  notation,  let  e.(dp  be  the  expected  earliness  given 
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di,  that  is,  c^d^  =  [  "x  fi(x  -  d.)  dx  , 

J  —CD 

and  Tj(d^)  be  the  expected  delay  given  d^,  that  is, 

Ti(di)  =  J0  X  fi*x  "  dP  **  * 

Then  the  combined  problem,  relaxing  scheduling  constraints  can  be 
formulated  as  follows: 

minimize  £  3HiTi  +  Aj/  T.  +  ^irici(di)? 

subject  to  £  $s./  T.  +  +  T.Cd^  /  T.  \  <  1  . 

Assuming  the  constraint  is  binding,  the  Lagrangian  relaxation  of  this 
problem  is 

t  -  l  +  A./  T.  +  V-e^d.)  +  e[5i/  T.  +  p.  +  r.(d.)/  T4]? 

where  0  £  0,  the  Lagrangian  multiplier  of  the  constraint,  measures  the 
imputed  value  of  machine  capacity.  The  optimal  lot  sizes  are 
determined  by  solving 


fr  -  Hi  -  rs  fAi +  9si +  eTi(di,]  -  0  • 

1  Ti 


Since  ■■  =  — s -  [A.  +  ©s .  +  0r.(d.)]  £  0  ,  L  is  convex  in  T. 


d2T.  T2  1 


i '  r 


Furthermore, 


Aj  +  0si  +  ©r^d.) 


Observe  that  this  is  similar  to  the  traditional  ELSP  solution  with  the 


expected  delay  being  added  to  the  setup  time.  The  optimal  delivery 
times  are  determined  by 


d  X 

d  d. 
1 


-hir.Fi(-d.)  +  §-  [1  -  F.(-d.)]  =  0 


d^  r  —  9 

Since  - 5  =  h.r.f.(-d.)  +  =-  f.(-d.)  £  0  ,  Z  is  convex  in  d.. 

d .  l 


Hence,  we  choose  di  such  that  F^(-d^)  = 


9 


9  +  h.r  *T. 

ill 


Observe  that  this  is  similar  to  the  solution  developed  in  previous 
sections  with  9  being  the  value  of  machine  capacity.  To  check 
convexity  in  both  variables. 


_si*_  .  -  [i  -  p  <-a  )]  . 

d  d.  d  T.  T.  1 

l  l  l 

The  determinant  of  the  Hessian  is 

fi  +  9si  +  +  f  ri(-diJ 

1  i  1 

1 


which  is  positive  if 


[A./9  +  s.  +  T.Cd^]  [1  +  h.r^/9]  >  1/2  [1  -  Fi(-di)32/  f i(-di)  . 


TV’ 


VT  ™ 


Observe  that  the  right  hand  side  of  this  inequality  gets  very  small  in 


the  right  tails  of  most  distributions.  For  example,  in  the  normal 


distribution,  at  the  75th  percentile  this  expression  equals  0.10,  at 


the  90th  percentile  it  is  0.03,  at  the  95th  percentile  it  is  0.012,  and 


at  the  99th  percentile  it  is  0.004  .  Since  we  are  focusing  on  the 


bottleneck  machine,  we  can  assume  that  deliveries  will  be  scheduled  to 


arrive  close  to  on-time,  hence  we  will  be  looking  at  the  right  tails  of 


the  distributions. 


4.5  A  Dynamic  Programming  Model. 


4.5.1  General  Formulation. 


The  traditional  ELSP  is  infeasible  if  the  capacity  of  the  machine 


is  exceeded,  that  is,  when  J  p.  >  1.  In  this  case,  it  is  impossible 
to  satisfy  the  total  demand.  Suppose,  however,  that  demand  is  dynamic 


rather  than  constant  and  that,  on  the  average,  demand  can  be  met. 


Consider  the  following  model  where  superscript  notation  is  used  to 


denote  specific  periods. 


Let  b^  be  the  backlog,  or  excess  capacity  required,  at  the  end  of 


period  j  —  normalized  by  machine  capacity,  that  is,  (hours  of 


backlog)/(hours  in  period).  We  assume  that  backlog  from  a  given  period 


uses  capacity  from  the  next  period.  Then  we  get  the  transition  equation 


J+1  -  J  = 


=  l  [p-  +  84/  +  T.(dJ+i)/  T.]  -  1 


bj+1  =  bj 


+  l  [p-+1  +  8./  t|+1  +  T.Cd^1)/  T.]  -  1 


T WJT 


'** 

i 

5 

J'i 

6 


« 

£ 


3 


Si 


Let  k  be  the  unit  cost  of  exceeding  machine  capacity.  The  penalty 
cost  due  to  backlog  (charged  at  the  end  of  each  period)  is 


,  '  ,  . 

P(bJ)  =  <  k  bJ  ,  for  bJ  >  0  , 

•  0  ,  for  bJ  £  0  . 

The  average  cost  for  period  j  is 

Cj  =  l  [1/2  h.rj(l  -  p|)Tj  4  h.rjc.Cdj)  +  K./  Tj]  . 


Let  GJ(bj_1-  bJ)  =  min  CJ 


subject  to  £  [p|  +  s./  t|  4  T.(dj)/  t|]  £  1  -  b^  14  b^  , 


where  b 


J"1  -  J  i 


bJ  is  the  relative  capacity  in  period  j  used  to  diminish 


the  backlog.  Then  we  can  formulate  the  problem  as  a  dynamic 
programming  problem. 

Let  (b^  *)  be  the  optimal  cost  of  going  to  the  end  of  the 
planning  horizon,  given  that  we  start  period  j  with  backlog  bJ  *.  For 
the  last  period,  assume  bn  =  0,  that  is,  we  have  no  backlog  at  the  end 
of  the  planning  horizon.  Then  for  the  last  period,  we  have 


Wn(bn_1)  =  Gn(bn  *) 


gm  11 

=  mm  C 

subject  to  l  [p"  4s./  T"  4  T.(dJ)/  T?  ]  £  1  - 

The  solution  to  this  is  given  by 


K.  4  0  l  [s.4  r.(dj)] 

1/2  h4rj  (1  -  Pp 


I  VO  «  1 ••J'.  ^  ' 


where  0  is  chosen  such  that  the  constraint  is  satisfied.  Observe  that 


larger  values  of  bn  *  make  the  constraint  tighter,  hence  the  imputed 
value  of  machine  time  0  becomes  higher.  This  causes  longer  cycle  times 
and  earlier  deliveries  and  thus  increases  Cn.  For  the  next  to  last 
period,  we  get 

wn-i(bn-2)  =  M.n  jP(bn-l)  +  Gn-1(bn_2-  bn-l^  +  ^^n-1^  ^ 
b""1 

We  only  need  to  consider  bn  ^  1  -  Z  p”  since  we  assume  bn  =  0.  In 
general,  the  following  recursion  exists: 

Wj(bj_1)  =  min|P(bj)  +  Gj(bj_1-  bj)  +  Wj+1(bj)$  . 
bj 

We  assume  that  the  boundary  conditions  are  bn  *  0  and  b^  =  0. 

4.5.2  A  Special  Case. 

Consider  the  following  special  case  of  the  general  problem 
developed  in  the  previous  section. 

Let  Ik  =  0  for  all  i,  that  is,  ignore  the  delivery  problem, 
and  let  K.  =  0  for  all  i,  that  is,  setup  costs  equal  zero. 

Then  for  the  solution  to  GJ(bJ-*~  b^)  we  get 


where  0  is  chosen  to  make  the  constraint  binding.  By  substitution  we 
get 


0  = 


i  si  j  »; 


JT^  -  jpj-  fj-1  *  »J] 


and 


Let  Bj  = 


I  *i  J 


* 

T-!  =  |  _ T _ . _ 

1  flT"  [1  -  l  pj  -  bJ_1  +  bJ] 


jJ  ^  _  b^  .  Then  we  can  write  G^  as 


i — l  . 


Gj(Bj)  = 


t  Z  J~h|~  ]  [  £  .j  J~hJ~  ] 

[»  -  z  o{  - 


GJ(BJ)  is  a  strictly  increasing  function  of  BJ  up  to  BJ  =  1  -  £  pj  and 
d  GJ/  d  BJ  =  GJ(BJ)  /  [1  -  Z  Pj  “  bJ]  • 

Since  BJ  is  restricted  to  be  less  than  1  -  ][  pj  ,  to  minimize  GJ 

choose  BJ  as  small  as  possible,  or  for  any  given  b^  ,  choose  bJ  as 
large  as  possible.  Looking  at  the  overall  problem  for  the  last  period, 
we  have 


Wn(bn  l)  =  Gj(bn_1) 


[  i  R~]  t  z  »t 

[1  -  X  p”  -  b"-1] 


For  the  next  to  last  period,  we  get 


Assuming  b 


2  0,  we  want  to  minimize 


n  1  n  1  n  1  C  Z  1  ]  [  £  S.  J  H°  1  ] 

Mn-1(bn  *)  =  k  bn  1  +  - - ^77 — 

[1  -  £  f>\  ~  bn'2  +  bn  *] 


[  ifilf]  [  l  *i  faf] 


[i  -  l  p”  -  b"'1] 


Observe  that  if  bn  *  <  0,  the  term  k  bn  *  drops  out. 


d  M 


n-1 


d  b 


n-1 


=  k  - 


C  £  J  Hf 1  ]  [  £  s.  J  H Y*  ] 


.n-1 


[1  -  £  p""1  -  bn~2  +  b11'1]  2 


c  i  [  i  fsf] 


[1  -  x  P*  -  b"'1]  2 


d2Mn_1 


2[  £  J  H”"1  ]  [  l  8i  j  Hf 1  ] 


d  b 


4 

n-1 


[i  -  i  p"-1  -  b”-2  +  b”-1]  3 


rT] 


tl  -  Z  p”  -  b"’1]  3 


2  o  . 


Hence  Mn  *  is  convex  in  b11  *  and  a  minimum  with  respect  to  bn  *  can  be 


n-1 


found.  Observe  that  if  b  £  0,  the  term  k  drops  out  of  the  first 


derivative  and  the  second  derivative  remains  the  same,  hence  convexity 
is  preserved. 
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Conditions  for  backlogging. 

Using  the  above  results,  we  can  then  define  the  condition  for 
backlogging  as 


k  < 


c  z  J 


h"'1  3  t  z  *i  J  h?'1  3 


[i  -  l  of1  -  b"-s  ]  2 


[  E  J~h?~3 

a  -  i  p"  ]  2 


For  example,  if  the  penalty  for  backlogging  k  were  very  small,  then  it 
might  be  advantageous  to  backlog  capacity  from  one  period  to  the  next. 
The  condition  for  negative  backlogging,  that  is,  getting  ahead  of 
demand,  is 


[  l  J 


H-_1  3  t  l  «i  J  H f 1 


[1  -  l  P 


n-1 


-  b"-2  ] 


This  can  be  defined  roughly  as  the  ratio  of  scaled  holding  costs  to 
slack  machine  capacity.  If  this  ratio  is  smaller  for  one  period  than 
for  the  subsequent  period,  then  it  is  advantageous  to  get  ahead.  The 
conditions  where  no  backlogging  is  optimal  can  be  defined  as  anything 
in  between  the  two  ranges  described  above.  Suppose  the  optimal  decision 
in  period  n-1  is  bn  *.  Then  for  period  n-2,  we  have 


Wn-2(bn  3)  -  min  jP(bn  2)  +  Cn~2(b”'3-bn"2)  +  2)j 

bn~2 


,  ,  n-a  .  I  [Zf-TJ  t£»t  3 

[i  -  i  P;-2  -  b"-3  +  b"-2] 
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9. 


F: 


k. 


C  l  J  n”'1  ]  [  l  »i  J  af 1  ]  [  l  J  ]  [  l  a,  J  H"  ] 


,n-l 


Up  n-1  .  n-2  .  n-ln 

-  E  p.  -  b  +  n  J 


[1  -  l  p”  -  b”'1] 


plus  an  additional  term  k  bn  1 ,  if  bn  1  >  0.  Observe  that  we  can 


define  ranges  for  backlogging  similar  to  that  in  period  n-1. 


4.6  Assembly  Model 


4.6.1  Introduction 

The  traditional  ELSP  assumes  constant,  continuous  demand.  An 
extension  to  this  is  suggested  by  the  OPT  problem  (5)  in  which  several 
parts  are  processed  through  a  machine  system  and  assembled.  In  this 
case,  inventory  holding  costs  for  IIP  are  charged  until  the  finished 
parts  are  assembled.  There  are  three  issues  that  must  be  addressed  in 
this  assembly  model: 


lot  sizing  (and  scheduling), 
sequencing, 

and  inventory  pattern. 


The  inventory  pattern  issue  will  be  discussed  in  the  following  two 
sections,  with  subsequent  comments  on  lot  sizing  and  sequencing. 


4.6.2  Two  Product  Case 

In  the  two  product  case  (see  Figure  4  -  6),  we  get  an  inventory 
pattern  as  illustrated  in  Figure  4-7. 


6i 


Assembly 
FIGURE  4-6 
Machine  Network 


FIGURE  4-7 

Two  Product  Inventory  Pattern 

In  this  inventory  pattern,  is  the  fraction  of  the  production  time  of 
i  during  which  inventory  is  accumulated.  During  the  other  portion 
(1  -  ap  of  the  production  time  all  parts  produced  are  used  for 
assembly. 

Note  that  pb(l  -  ab)pfaT  =  Pa«aPaT  and  P&(1  -  «a>PaT  =  Pb«bPbT  • 

The  cycle  length  is  defined  by 

T  =  P,T  +  PbT  +  Sa  +  Sb 

■  <S,  +  V  /  <»  -  P,  -  0b>  • 

Si  2  s.  +  Tj(d^)  ,  for  i  =  a  and  b. 


where 


Then  we  can  formulate  an  expression  for  the  average  inventory,  St,  of 
each  part  as  follows: 


SA  =  (paVaT  /  2)(VaT  1  T)  +  (paaapaT)(Sb  1  T) 

+  (pa“aPaT  1  "  “a^b7  1  T) 

=  (Pa  /  Pb)  «aPa  /  2  [VaT  (pa+  pb)  +  2  pbSb] 

*8  =  (pbVbT  /  2)(abpbT  7  T)  +  (pbVbT)(Sa  1  T) 

+  (PbKbPbT  /  2)[d  -  «a)paT  /  T] 

=  (Pa  /  Pb)(l  -  «ft)pa/2  1(1  -  «a)paT  (pa+  pb)  +  2  PbSj  . 
The  total  average  cost  to  minimize  is 

C  =  Ha(pa  7  Pb)  Va  1  2  CVaT  (pa+  pb)  +  2  pbSb] 

+  Hb(pa  '  pb)(l  "  aa)pa  1  2  [(1  "  aa)paT  (pa+  pb)  +  2  pbSa] 


and  similarly  for  B. 


jl 


.%  a> 


This  gives  T  =  !>a(da)  +  sa  +  Tb(db>  4  sj  /  [1  -  p.  -  pj  - 


a  KbJ 


Now  that  we  know  how  to  choose  S  and  S,  ,  we  can  find  the  minimum 

a  b 

with  respect  to  «a« 

d  C  /  d  «a  =  Ha(pa  /  pb)  Pa  /  2  [«apaT  (pa4  pj  +  2  p^Sj 


a  v*a  *V 


+  Ha(pa  /  pb)  «apa  /  2  [paT  (pa+  pb)] 


-  Hb(pa  /  pb)  Pa  /  2  [(1  -  «JpJ  (p^+  pK)  +  2  pKSj 


a^a‘  Via  *V 


-  Hb(Pa/  Pb)d  -  OP,  /  2  pj  (p.4  pj 


a  Via 


=  papa  /  pb  ^Va7  (pa+  pb>  +  pbSb^ 

-  Hb[(l  -  «a)paT  (pa+  pb)  4  pbSa]? 

d2C  /  d  ora2  =  Pa Pa/Pb  (Ha4  Hb)(pa4  pb)  paT  £  0 

Hence  C  is  convex  in  a  .  The  optimal  value  of  a  is 

a  a 

*  Hh  Pu  H.  S  H  S. 

«  =  _ L_  +  b  5  -  b  _1_  _  _ L_  JL  l 

a  H  4  H,  p  4  p.  !  H  4  H,  p  T  H  4  H,  p  T 
a  b  a  *d  a  b  *a  a  b  Ka 

Suppose  that  H  =  H,  and  S  =  S  .  Then  a  =  1/2.  This  gives  us  a 
general  solution  to  the  two  product  case.  However,  an  important 
question  is  when  does  one  product  dominate  the  problem  so  that  we  never 
carry  any  inventory  of  that  product? 

Conditions  for  dominant  product. 


Under  what  conditions  should  we  never  accumulate  finished  inventory  of 
* 

A,  that  is,  a  £  0?  These  conditions  can  be  summarized  as 

a 


.-vivo, 
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ja  ,  V  *b  +  V  T 
Hb  sb  7  T 


P  +  Pv 


1  -  p  -  p, 

Ka  *d 


This  condition  holds  when  the  holding  cost  of  A  is  much  larger  than  the 
holding  cost  of  B.  The  last  inequality  is  strict  unless  S&  =  0.  For 
example,  if  the  machine  utilization,  p&  +  p^,  equals  0.8,  then  H&/Hb 
must  be  at  least  4.  To  find  the  minimum  with  respect  to  d  , 

d  C  /  d  da  =  Hbpa(l  -  «a)pa  [1  -  Fa(-da)]  -  r  haFa(-da)  -  0. 

dS  C  /  d  da2  »  Hbpa(l  -  «a)Para(-da)  +  r  hafa(-da)  i  0 


Hence  C  is  convex  in  d  .  Choose  d^  such  that 

a  a 

H  p  (1  -  «  )p 

F  (-d  )  =  - — - — - 

a  a  Hp,(l  -  «  >P*  +  r  K 


a‘  a 


a  r  a 


To  check  convexity  in  d  and  «  , 

a  a 


d  0  /  d  da  d  da  =  -Hbpapa[l  -  Fa(-da)]  . 
The  determinant  of  the  Hessian  is  positive  if 


P.P,  /  »b  <V  Hb)(',,+  "b>  paT  [V.(‘  -  %>p»fa(-d»)  +  r  Va<-da)] 

Z 


*  HbV  -  Fa(-da>] 


or 


r  h  [1  -  F  <— d  )  J 

»b  P»TS  *  (1  '  %)  +  Vapa  *  "  V-da> 


H  +  H,  p  +  p. 

a  b  ra  rb 


Observe  that  the  right  hand  side  o^  his  inequality  is  the  same  as  that 
developed  in  Section  4.4,  and  convexity  again  holds  for  reasonable 
problems. 


$ 


Three  Product  Case 


Consider  three  products  produced  on  a  rotation  cycle  (see 
Figure  4-8). 


FIGURE  4-8 

Three  Product  Inventory  Pattern 


The  average  inventory  of  each  part  is 

&A  =  <paVaT  /  2>  <VaT  '  T)  +  (paVaT)  <V  T) 

+  [paVaT  1  2  +  pc(1  ~  ac)pcT  1  2]  [(1  '  ab)pbT  /  T] 
+  [P„(l  -  «  )p  T]  [(«.  p.T  +  SJ  /  T] 


+  [p„(l  "  «  )P  T  /  2]  [(1  -  Op  T  /  T] 


*b  =  [pb(i  -  «a)p  t  /a2]a[(i  -  OPT  /aT]a  +  (PbvbT)  <sc  '  T> 

+  £<p  «  PbVb2>  +  p  <1  I  «  >p  l 

•  [(1  -  ac)pcT  /  T]pa(l  -  «a)paT]  [(«CPCT  +  Sc)  /  T] 

Sc  =  [(PC«CPCT  /  2)  +  pb(l  -  «b)pbT  /  2]  [(1  -  «a)paT  /  T] 

+  [pb( 1  -  «b)pbT  /  2]  [(«bPbT  +  Sb)  /  T] 

+  [Pb(l  -  «b)pbT  /  2]  [(1  -  «b)pbT  /  T] 

+  (pcacpcT  /  2)  (acpcT  /  T)  * 

The  following  is  the  total  inventory  cost: 

C  =  H  a  +  H.  SL  +  H  S  . 

a  a  b  d  c  c 


C  is  convex  in  ot^  and  the  determinant  of  the  Hessian  with  respect  to  a. 

and  a.  is  positive,  hence  we  have  convexity  in  two  variables. 

J 


P  +  Pv 

*a 


(1  -  of  )p  + 


C/rC  n  +  n,  n  T 
ra  ^b  Ka 


,  -H  S,  +  2H,S  -HS, 

1  <  a  b  be  c  b  ? 

«'pc  U  _L  U  -J-  U  > 


H  +  H,  +  H 
a  b  c 


Ha  +  H.  +  H 
a  b  c 


u  + 


pa  pb 


VVV 

Suppose  pa=  pb  =pc  ,  Sa=  Sb=  Sc  ,  and  Ha=  Hb=  Hc>  Then  ora  =  2/3  . 


4.7.4  Sequencing  When  One  Product  is  Dominant. 

Suppose  one  product,  call  it  product  A,  is  dominant  so  that  no 


inventory  of  that  product  is  ever  accumulated.  How  should  the  remaining 
products  be  sequenced?  Consider  a  three  product  case  with  inventory 
pattern  shown  in  figure  4-9. 


2- 


Product  A 
Inventory 


FIGURE  4-9 

Three  Product  Inventory  Pattern 

The  average  inventory  for  B  is 

%  =  (rT  /  2)  (pbT  /  T)  +  r  (Sg  /  T)  +  (rT  /  2)  (p#T  /  T)  . 

The  average  inventory  for  C  is 

Sc  -  (rT  /  2)  (pcT  /  T)  +  r  (Sa  /  T)  +  (rT  /  2)  (p&T  /  T) 

+  r  [(Sb  +  pbT)  /  T]  . 

Observe  that  the  only  term  in  these  expressions  that  is  dependent  on 

the  sequence  of  B  and  C  is  r  [(Sb  +  pbT)  /  T].  Then  using  the 
sequence  dependent  costs,  B  should  follow  C  if 

Hcr  [(Sb  +  pbT)  /  T3  *  Hbr  [<»c  +  pcT)  /  T] 
or  1/Hb  (Sfe  +  pbT)  ^  1/Hc  (Sc  +  pcT)  . 

This  is  similar  to  classic  scheduling  theory's  weighted  shortest 
processing  time  rule  except  that  here  we  use  weighted  longest 

processing  time  with  the  weights  being  the  inverse  of  holding  costs. 


CHAPTER  5 


AGGREGATE  PLANNING  ALGORITHM 

5.1  Algorithm  Development 

5.1.1  Background 

Section  4.4  presented  a  model  which  combined  the  issues  of 
deliveries  to  a  bottleneck  machine  and  lot  sizing.  The  optimal  order 
intervals  and  requested  delivery  times  were  as  follows: 

Ai+e(Si+T.(d.)) 

Ti 

and 

F,(-di)  *  -2- -  ,  (1) 

9+h.r^ 

where  8,  the  imputed  value  of  machine  capacity,  is  chosen  such  that  the 
machine  capacity  constraint  is  satisfied  at  equality,  that  is, 

I$S./T.  4  p.  4  T.(d.)/T.S  =  1 
or 

l\Si  4  T.(d.)$/T.  =  1  -  Spi.  (2) 

If  the  setup  cost  is  zero,  that  is,  the  value  of  setup  consists  only  of 
the  value  of  lost  machine  capacity. 
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S.+T.(d.) 

Ti-r«  -4^- 


Substituting  this  into  the  capacity  constraint  (2),  ire  have 


s  1  WV 


=  1  -  2p. 


Z  H. (S.+t. (d.J 
i v  i  i v  i; 

l-zp. 


This  leads  to  the  following  algorithm  which  iterates  on  the 
unknown  imputed  value  of  machine  capacity: 


Step  1.  Let  T^(d^)  =  0  for  all  i. 


Step  2.  Compute  6  using  (4). 


z  H.(s.+T.(d.; 
x x  x  i v  i ; 


Step  3.  Compute  for  all  i  using  (3). 
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Step  4.  Compute  for  all  i  using  (1). 


F.(-d.)  =  - 

1  i  «,r 


e 


"Wi 

The  inverse  of  F(.)  depends  on  the  form  of  the  distribution 
distribution  of  delivery  tiroes. 


Step  5.  Compute  T^(d^)  for  all  i.  t  (.)  depends  on  the  form  of 
the  distribution  of  delivery  times. 

Step  6.  Go  to  step  2  until  0  converges.  Convergence  is  assured  by 
the  convexity  verified  in  chapter  4. 


5.1.3  Normal  Distribution  Assumption 

If  we  assume  that  the  delivery  times  to  the  bottleneck  machine  are 
normally  distributed,  for  step  5  of  the  algorithm  we  get 

Ti(di)  =  o.  f(-d  ./op  +  d.  F(di  /a.) 

where  f(.)  is  the  standard  normal  density  function  and  F(.)  is  the 
standard  normal  cumulative  distribution  function.  This  algorithm  is 
coded  in  PASCAL  in  Appendix  1. 
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K 


5.2  Results 

Two  key  questions  to  ask  using  this  model  are: 

When  can  we  ignore  the  delivery  issue  in  scheduling  a  machine  and 
what  affect  do  the  parameters  of  this  model  have  on  the  overall 
system?  These  questions  can  be  answered  by  looking  at  the  results  of 
various  data  input  to  the  model. 


5.2.1  Holding  Costs  for  Raw  Parts 

We  can  assume  that  the  ratio  of  holding  costs  for  raw  parts  to  the 
holding  costs  for  processed  parts  is  between  0  and  1,  that  is,  the 
value  of  the  raw  parts  lies  somewhere  between  zero  and  the  value  of  the 
processed  parts.  We  can  compare  the  resulting  lot  sizes  determined  from 
the  aggregate  model  with  the  corresponding  lot  sizes  determined 
ignoring  the  delivery  issue.  The  ratio  of  these  lot  sizes,  call  it  the 
.3-ratio,  is  one  when  the  value  of  raw  parts  is  zero,  hence  we  can 
ignore  the  delivery  issue  in  scheduling  a  machine  when  the  relative 


Pi 


I? 


t? 


value  of  raw  parts  is  small.  In  this  case,  we  can  carry  sufficient 
inventory  of  raw  parts  to  keep  the  bottleneck  machine  fully  utilized 
because  the  raw  parts  inventory  is  very  cheap.  However,  Figure  5-1 
shows  that  as  the  relative  value  of  raw  parts  increases,  the 
corresponding  lot  sizes  in  the  aggregate  model  also  increase.  Hence,  as 
the  relative  value  of  raw  parts  increases,  it  becomes  more  important  to 
consider  the  delivery  issue  when  scheduling  a  machine. 


Part 


Data  Used  in  Figure  5-1 


1 

2 

3 


pi 

0.25 

0.33 

0.33 


8i 

3 

2 

1 


l 

60 

40 

20 


h. 

l 

4 

2 

3 


FIGURE  5-1 
Graph 


h/h 


5.2.2  Variance  of  Delivery  Times 

If  the  variance  of  delivery  times  is  small,  we  can  ignore  the 
delivery  issue  when  scheduling  a  machine.  To  quantify  how  small  the 
variance  should  be,  we  can  compare  the  standard  deviation  of  delivery 
time  to  the  setup  times.  As  shown  in  Figure  5-2,  if  the  standard 
deviation  of  delivery  time  is  small  relative  to  the  setup  times,  we  can 


ignore  the  delivery  issue.  We  also  notice  from  Figure  5-3  that  the 
total  system  costs  go  up  at  an  almost  linear  rate  with  respect  to  the 
standard  deviation  of  delivery  time.  Hence  the  variance  of  the  input 
process  plays  a  critical  role  in  both  lot  sizes  and  overall  costs. 

TABLE  5-2 


FIGURE  5-2 
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Graph 

Reducing  Setup  Times 

Much  of  the  literature  purporting  the  advantages  of  the  Japanese 
philosophy  of  scheduling,  just-in-time,  support  reducing  setup  times  on 
machines  as  a  means  to  decrease  lot  sizes.  However,  if  we  look  at 
reducing  setup  times  in  the  aggregate  model,  we  find  that  the  variance 
of  the  input  process  plays  a  critical  role  in  determining  how  much  the 
lot  sizes  can  be  reduced.  As  shown  in  Figure  5-4,  reducing  setup  times 
produces  a  corresponding  reduction  in  lot  sizes  until  the  setup  time 
gets  small  compared  to  the  variance  of  the  input  process.  In  Figure  5- 
5,  we  see  that  the  ratio  of  lot  sizes  determined  from  the  aggregate 


model  compared  to  the  corresponding  lot  sizes  determined  ignoring  the 
delivery  issue  increases  almost  exponentially  as  the  setup  times 
approach  zero.  Hence,  in  trying  to  achieve  just-in-time  via  small  lot 
sizes,  the  variance  of  the  input  process  must  be  reduced  as  well  as 
reducing  the  setup  time  on  the  machine. 

TABLE  5-3 

Data  Used  in  Figures  5-4  and  5-5 


Part 

pi 

s . 

X 

a. 

l 

h. 

l 

1 

0.25 

3 

60 

4 

2 

0.33 

2 

40 

2 

3 

0.33 

1 

20 

3 

FIGURE  5-4 


CHAPTER  6 


CONCLUSIONS 

If  we  put  the  lot  scheduling  problem  in  the  context  of  the 
hierarchy  of  decision  models  and  look  both  at  the  physical  setting  and 
the  nature  of  the  problem,  we  can  get  better  solutions  to  a  given 
problem  in  a  more  realistic  setting.  Due  to  the  inherent  difficulty  in 
solving  scheduling  problems,  the  best  we  can  usually  hope  for  is 
heuristics  which  provide  good  solutions  over  a  reasonable  range  of 
problems.  However,  if  we  encounter  a  problem  with  a  dominant  product, 
we  can  make  use  of  that  fact  to  simplify  the  search  for  a  solution. 

The  concept  of  a  dominant  product  can  be  used  to  determine  when 
to  insert  idle  time  into  a  schedule.  In  this  case,  a  dominant  product 
is  one  with  dominant  holding  costs.  If  in  a  given  problem  we  have  a 
dominant  product  and  the  remaining  products  have  low  machine 
utilization,  that  is,  we  have  slack  machine  time  available,  then  we 
can  produce  a  better  schedule  by  inserting  idle  time  for  the  dominant 
product. 

The  role  of  the  dominant  product  also  tells  us  when  not  to  use 
the  zero  switch  rule  (ZSR).  We  can  conclude  that  the  zero  switch  rule 
(ZSR)  is  a  good  scheduling  policy  for  most  problems.  The  exception  to 
this  happens  when  the  ZSR  solution  yields  lumpy  production  patterns 
for  a  dominant  product.  In  this  situation,  we  can  sometimes  improve 
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upon  the  ZSR  solution  by  making  the  dominant  product's  production 
pattern  more  even  while  using  a  non-zero  switch  on  another  product. 

In  looking  at  the  context  of  a  machine  being  scheduled,  if  there 
is  variability  in  the  input  process  to  the  machine,  we  can  achieve 
lower  overall  system  costs  by  considering  the  issues  of  determining 
delivery  times  and  lot  sizes  concurrently  instead  of  looking  at  each 
independently,  lfe  can  use  this  approach  both  in  aggregate  planning  and 
in  real-time  detailed  scheduling.  When  determining  the  optimal 
delivery  times,  we  balance  the  cost  of  holding  raw  parts  with  the 
value  of  lost  machine  capacity  in  a  manner  similar  to  the  Newsboy 
Problem.  This  approach  can  be  extended  to  consider  problems  with  non¬ 
constant  demand  when  we  allow  backlogging  of  demand  and  machine 
capacity.  However,  we  find  that  much  of  the  time,  we  will  still  try  to 
satisfy  current  period  demand  with  current  period  production. 

The  results  of  the  aggregate  model  combining  the  issues  of 
delivery  and  lot  sizes  show  that  the  conditions  when  we  should 
consider  the  delivery  issue  in  conjunction  with  the  lot  size  issue  are 
when  the  holding  cost  of  raw  parts  is  high  with  respect  to  the  holding 
cost  of  processed  parts  and  when  the  variance  of  delivery  times  is 
high  with  respect  to  the  corresponding  setup  times.  We  find  that  the 
overall  system  costs  increase  at  an  almost  linear  rate  with  respect  to 
the  standard  deviation  of  delivery  times. 

In  summary,  we  find  that  by  placing  a  given  problem  in  its  proper 
context,  we  can  more  effectively  derive  solutions  to  the  real  problem 


at  hand. 


APPENDIX  1 


AGGREGATE  PLANNING  ALGORITHM 


PROGRAM  sp  (INPUT, OUTPUT, mdata);  $***  TABLE  OF  CONTENTS 

1  Declarations 

2  Readin  —  reads  input  data 

3  Setup  —  sets  up  in  proper  format 

4  Norminv  —  finds  inverse  of  normal  CDF 

5  Normden  —  calculates  density  of  std  normal  rv 

6  Cum  —  computes  value  of  the  normal  CDF 

7  Thetacalc  —  calculates  value  of  machine  time 

8  Results  —  summarizes  part  schedules 

9  Display  —  writes  results 

10  CONTROLLING  PROGRAM  .«***»».«**»,,****** 

CONST  p  =  3;  5  number  of  parts  £ 

n  =  4;  ^number  of  iterations! 

lead  =  5;  \ lead  time  for  delivery! 


TYPE  index  =  integer; 

parray  =  array  [l..p]  of  real; 


VAR  mdata;  text;  \ input  file! 

r  :  parray;  ^demand  rate! 

su  :  parray;  \ setup  time! 

pr  ;  parray;  ^production  rate! 

ru  ;  parray;  ^relative  utilization  -  r/pr! 

h  :  parray;  ^holding  cost  after  processing! 

hb  :  parray;  ^holding  cost  before  processing! 

sd  :  parray;  \ standard  deviation  of  delivery  time! 

t  :  parray;  \ cycle  length  -  time  between  production  runs 

fcum  :  real;  ^fraction  of  deliveries  on  time! 

comp  :  real;  jl  -  fcum! 

ninv  :  real;  formal  inverse  of  comp! 

fden  :  real;  ^density  of  delivery  dist'n! 

d  :  parray;  ^scheduled  delivery,  neg  =  early,  pos  =  late! 

tau  :  parray;  ^expected  delivery  delay! 

mtau  ;  parray;  ^expected  early  delivery! 

suml  :  real;  )sum  of  sqrt(h)  *  (su  +  tau)  ! 

sum2  :  real;  jsum  of  su  +  tau  ! 

si  :  real;  jslack  =  1  -  sum  of  ru  ! 

theta:  real;  jvalue  of  machine  time  =  suml**2/(sum2*l**2) ! 
cost  :  real;  )cost  of  solution  at  last  iteration! 
i,j  :  index; 


PROCEDURE  readin;  Pread  in  data? 

var  i:  index; 
begin  Preadin? 
reset(mdata) ; 
for  i  :=  1  to  p  do  begin 
read(mdata,r[i]) ; 
read(mdata,su[i]) ; 
read(mdata,pr[i]); 
read(mdata,sd[i]) ; 
read(mdata,h[i]) ; 
read(mdata,hb[i]) ; 
end; 

close(mdata) ; 
end;  preadin? 


PROCEDURE  setup;  p setup  data  in  proper  format ? 

var  i:  index; 
begin  | setup ? 
si  :=  1.0; 

for  i  :=  1  to  p  do  begin 
ru[i]  :=  r[i]  /  pr[i]; 
si  :=  si  -  ru[i]; 

tau[i]  :=  0;  P intitializes  tau  to  zero? 
end; 

writelnC  slack  capacity  =  ',sl); 

end;  P setup? 

PROCEDURE  norminv; \ computes  the  appoximate  inverse  of  the  normal  CDF? 

preference:  HANDBOOK  OF  MATHEMATICAL  FUNCTIONS  by  ? 


PAbramowitz  and  Stegun,  error  <  .00045  ? 
P given  comp,  where  0  <  comp  <  .5  ,  finds  ? 
P complementary  cumulative  inverse  ninv,  ? 
P where  1  -  F(x)  *  comp  ? 


var  tl,t2,nl,dl,x:  real; 
begin  p norminv? 
tl  :=  1.0; 

if  comp  >  0.5  then  begin 
tl  :=  -1.0; 
comp  :=  1.0  -  comp; 
end; 

t2  :=  sqrt(2*ln(l. 0/comp)) ; 
nl  :=  2.515517  +  0.802853U2  +  0.010328*t2*t2; 
dl  :=  1.0  +  1.432788*t2  +  0. 189269*t2*t2  +  0.001308*t2*t2*t2; 
x  :=  t2  -  nl/dl; 
ninv  :=  tl*x; 
end;  P norminv? 
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PROCEDURE  normden;  $  calculates  density  of  std  normal  rv  $ 

var  n2:  real; 
begin  \ normden ? 

n2  :=  -(ninv*ninv)/2.0; 
fden  :=  0.3989423  *  exp(n2); 
end;  \ normden ? 

PROCEDURE  cum;  P computes  the  appoximate  value  of  the  normal  CDF  ? 

preference:  HANDBOOK  OF  MATHEMATICAL  FUNCTIONS  by? 
lAbramowitz  and  Stegun,  error  <  .00001  ? 

P given  ninv,  finds  value  of  the  normal  CDF,  fcum  ? 
var  tl,t2:  real; 
begin  Pcum? 

if  ninv  >  0.0  then  begin 

tl  :=  1.0  /  (1.0  +  0.33267*ninv) ; 

t2  :=  0.4361836*tl  -  0.1201676*tl*tl  +  0.937298*tl*tl*tl; 

fcum  :=  1  -  fden*t2; 

end 

else  begin 

tl  :=  1.0  /  (1.0  -  0.33267*ninv) ; 

t2  :=  0.4361836*tl  -  0.1201676*tl*tl  +  0.937298*tl*tl’*tl; 
fcum  :=  fden*t2; 
end; 
end;  pcum? 


PROCEDURE  thetacalc;  Pcalculates  value  of  macine  capacity? 
var  i:  index; 
begin  Pthetacalc? 
suml  :=  0.0; 
sum2  :=  0.0; 

for  i  :=  1  to  p  do  begin 

suml  :=  suml  +  sqrt(h[i]>  *  (su[i]  +  tau[i]); 

sum2  :=  sum2  +  su[i]  +  tau[i]; 

end; 

theta  :=  suml*suml  /  (sl*sl); 
end;  Pthetacalc? 


PROCEDURE  results;  Psummarize  results? 

var  i:  index; 

cut:  real;  ptruncated  left  tail? 

begin 

cost  :=  0.0; 

for  i  :*  1  to  p  do  begin 

if  lead  +  d[i]  -  3*sd[i]  <  0.0  then  begin 

ninv  :=  -lead  /  sd[i];  Ptruncates  delivery  to  lead  time? 
normden;  pcomputes  density  of  normal  rv? 


cum;  ^computes  cumulative  dist'n  of  ninv? 

cut  :=  sd[i]  *  fden  +  lead  *  fcum;  ^expected  early? 

mtau[i]  :=  mtau[i]  -  cut; 
end; 

cost  :=  cost  +  h[i]*t[i]  +  hb[i]*r[i]*mtau[i] ; 
end; 

end; 


PROCEDURE  display;  Jwrite  results? 

var  i:  index; 
begin  ^display? 

writelnC 'INPUT  DATA'); 

writelnC'  r/p  s  sd  h  hb'); 

for  i  :=  1  to  p  do  begin 

write ln(ru[i] ,su[i] ,sd[i],h[i] ,hb[i]) ; 
end; 

writeln( 'RESULTS'); 

writeln( 'Part  No  Order  Int  Delivery  Time  Expected  Delay'); 

for  i  :=  1  to  p  do  begin 

write lnCi,t[i] ,d[i] ,tau[i]) ; 
end; 

writelnC  COST  OF  SOLUTION  =  ’.cost); 
end;  ^display? 


begin  $sp? 

writelnC 'Start  program'); 

readin;  $read  in  data? 

setup;  \ setup  data  in  proper  format? 

writeln; 

for  i  :=  1  to  n  do  begin 

thetacalc;  ^calculate  value  of  machine  time? 

writelnC'Iteration',i, '  theta  =  ', theta); 

for  j  :=  1  to  p  do  begin 

t[j]  :=sqrtCtheta  *Csu[j]+tau[j])/h[j]) ;  ^production  interval? 
fcum  :=theta  /  Ctheta+r[j]*hb[j]*t[j]) ;  ^deliveries  on  time? 
comp  :=  1.0  -  fcum;  \ complementary  cumulative  of  fcum? 

norminv;  ^computes  functional  inverse  of  fcum? 

d[j]  :=  “  ninv  *  sd[j];  $  scheduled  delivery? 

normden;  J computes  density  of  normal  rv? 

tau[j]  :=  sd[j]  *  fden  +  d[j]  *  Cl  -  fcum);  ^expected  delay? 
mtau[j]  :=  sd[j]  *  fden  -  d[j]  *  fcum;  ^expected  early? 

end; 
end; 

results;  ^summarize  results? 

display;  jwrite  reslults? 

end.  fsp? 


APPENDIX  2 

THE  BAKER  AND  BOMBERGER  PROBLEMS 

The  Bomberger  problem  (3)  is  defined  by  the  following  parameters: 


i 

s . 

l 

Ai 

hi 

Pi 

pi 

Ti 

1 

.125 

15 

.00065 

30000 

400 

167.5 

2 

.125 

20 

. 01775** 

8000 

400 

37.7 

3 

.25 

30 

.01275 

9500 

800 

39.3 

4 

.125 

10 

.01 

7500 

1600 

19.5 

5 

.5 

110 

.2785 

2000 

80 

49.7 

6 

.25 

50 

.02675 

6000 

80 

106.6 

7 

1.0 

310 

.15 

2400 

24 

204.3 

8 

.5 

130 

.59 

1300 

340 

20.5 

9 

.75 

200 

.09 

2000 

340 

61.4 

10 

.125 

5 

.004 

15000 

400 

39.3 

Ti*  =  LAi/H±31/2 

"•Originally  .01175,  however  all  subsequent  authors  have  used  .01775. 


For  this  problem,  we  get  the  following  results: 


l  Pi  =  0.88 

i 

I  l  Si/T*  =  0.07 

[ 

* 

£  s./T.  +  p.  =  0.95  (reference  constraint  (2-2) 

1  in  Chapter  2). 


Hence,  only  95%  of  the  machine  capacity  is  utilized  in  the  theoretical 
optimum  solution. 
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The  Baker  problem  (1)  is  defined  by  the  fol lowing  parameters: 


8i 

A. 

l 

hi 

Pi 

pi 

Hi 

Ti 

.08 

75 

.01 

2500 

.08 

.92 

9.03 

.04 

30 

.10 

1000 

.25 

9.38 

1.79 

.02 

25 

.04 

500 

.20 

1.60 

3.95 

.12 

35 

.08 

200 

.35 

1.82 

4.39 

T.*  =  [A./H.]l/2 
l  L  r  iJ 

For  this  problem,  we  get  the  following  results: 
l  p.=  0.88 

l  s./T*  =  0.06 

s./T.  +  £  p.  =  0.94  <  1  (reference  constraint  (2-2) 

in  Chapter  2). 


Hence,  only  94%  of  the  machine  capacity  is  utilized  in  the  theoretical 
optimum  solution. 
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